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ABSTRACT 

We study the elemental abundances of C, N, O, Al, Si, S, Cr, Mn, Fe, Ni, and 
Zn in a sample of 14 damped Lya systems (galaxies) with H I column density A^(H 
I)> 10^'^ cm~^, using high quality spectra of quasars obtained with the 10m Keck 
telescope. To ensure accuracy, only weak, unsaturated absorption lines are used to 
derive ion column densities and elemental abundances. Combining these abundance 
measurements with similar measurements in the literature, we investigate the chemical 
evolution of damped Lya galaxies based on a sample of 23 systems in the redshift 
range 0.7 < z < 4.4. The main conclusions are as follows. 

1. The damped Lya galaxies have (Fe/H) in the range of 1/10 to 1/300 solar, clearly 
indicating that these are young galaxies in the early stages of chemical evolution. The 
N(Ji I)-weighted mean metallicity of the damped Lya galaxies between 2 < z < 3 is 
(Fe/H)=0.028 solar. There is a large scatter, about a factor of 30, in (Fe/H) at z < 3, 
which we argue probably results from the different formation histories of the absorbing 
galaxies or a mix of galaxy types. 

2. Comparisons of the distribution of (Fe/H) vs redshift for the sample of damped 
Lya galaxies with the similar relation for the Milky Way disk indicate that the damped 
Lya galaxies are much less metal-enriched than the Galactic disk in its past. Since 
there is evidence from our analyses that depletion of Fe by dust grains in the sample 
galaxies is relatively unimportant, the difference in the enrichment level between the 
sample of damped Lya galaxies and the Milky Way disk suggests that damped Lya 
galaxies are probably not high-redshift spiral disks in the traditional sense. Rather, 
they could represent a thick disk phase of galaxies, or more likely the spheroidal 
component of galaxies, or dwarf galaxies. 

3. The mean metallicity of the damped Lya galaxies is found to increase with 
decreasing redshift, as is expected. All four of the damped Lya galaxies at z > 3 in our 
sample have (Fe/H) around 1/100 solar or less. In comparison, a large fraction of the 
damped Lya galaxies at z < 3 have reached ten times higher metallicity. This suggests 
that the time around z = 3 may be the epoch of galaxy formation in the sense that 
galaxies are beginning to form the bulk of their stars. Several other lines of evidence 
appear to point to the same conclusion, including the evolution of the neutral baryon 
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content of damped Lya galaxies, the evolution in the quasar space density, and the 
morphology of z > 3 galaxies. 

4. The relative abundance patterns of the elements studied here clearly indicate 
that the bulk of heavy elements in these high-redshift galaxies were produced by Type 
II supernovae; there is little evidence for significant contributions from stellar mass loss 
of low-to-intermediate mass stars, or from Type la supernovae. 

5. Although earlier studies have attributed the overabundance of Zn relative to Cr 
or Fe found in damped Lya galaxies to selective depiction of Cr and Fc by dust grains, 
such an interpretation is inconsistent with many of the other elemental abundance 
ratios seen in these galaxies, most notably N/0 and Mn/Fe. Several other tests also 
indicate that there is no significant evidence for dust depletion in these galaxies. We 
suggest that the overabundance of Zn relative to Cr in damped Lya galaxies may be 
intrinsic to their stellar nucleosynthesis. If this interpretation is correct, it will provide 
important new information to the theory of stellar nucleosynthesis. 

6. The absorption profiles of Al III in damped Lya galaxies are found to resemble 
those of the low-ionization lines. The profiles of Si IV and C IV absorption, while 
resembling each other in general, are almost always different from those of the 
low-ionization absorption lines. These results suggest that Al III is probably produced 
in the same physical region as the low-ionization species in the absorbing galaxies, 
while the high-ionization species (Si IV and C IV) mostly likely come from distinct 
physical regions. 

7. We discuss possible ways to obtain information on the history of star formation 
(i.e., continuous or episodic) in damped Lya galaxies, and on the shapes of the stellar 
initial mass functions. 

8. We review the evidence for, and against, the hypothesis that damped Lya 
galaxies are disks or proto-disks at high rcdshifts, and discuss the implications. 

9. We determine upper limits on the temperature of the cosmic microwave 
background radiation at several redshifts using absorption from the fine structure level 
of the C II ion. These upper limits are consistent with the predicted increase of Tcmb 
with redshift. 

Subject headings: cosmology: cosmic microwave background - early universe - galaxies: 
abundances - galaxies: ISM - nucleosynthesis - quasars: absorption lines 

1. INTRODUCTION 

A fundamental question in the study of galaxy formation and evolution is: what is the 
chemical evolutionary history of galaxies? While some information pertaining to the general 
question of galactic chemical evolution has been obtained from studies of stellar populations and 
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the interstellar medium in the Milky Way and in nearby external galaxies, direct observations of 
the chemical evolution of galaxies in the early universe has been lacking. The situation changed 
when it was recognized that the narrow absorption lines in the spectra of high-redshift quasars 
can be used to probe conditions in the early universe. One of the most important developments 
along this direction has been the recognition of "damped Lya absorption systems" as a significant 
population of quasar absorbers and the subsequent systematic studies of their basic properties, 
which demonstrate that damped Lya galaxies are likely to be the progenitors of present-day 
galaxies (Wolfe et al. 1986; Wolfe 1995; Lanzetta, Wolfe, & Turnshek 1995). 

The pioneering work by Meyer, Welty, &; York (1989) and by Pettini, Boksenberg, & Hunstead 
(1990) clearly demonstrated that one can extract significant information on the abundances of 
chemical elements in damped Lya galaxies through the analysis of the absorption lines in quasar 
spectra. Pettini et al. (1994) conducted the first systematic survey of elemental abundances in 
damped Lya galaxies, though, for practical reasons, the study was limited to the elements Zn and 
Cr. Nonetheless, their study revealed the important fact that damped Lya galaxies are generally 
metal poor (mean (Zn/H)~ 1/10 solar), confirming that these are young galaxies in the early 
stages of chemical enrichment. It was also noted that the element of Cr is generally less abundant 
than Zn relative to their solar proportion. Since Zn and Cr are thought to be made in solar 
proportions throughout the Milky Way based on their abundances in Galactic disk and halo stars 
(see discussion in §4.3), the overabundance of Zn relative to Cr found in damped Lya galaxies was 
taken as evidence that some Cr was depleted from the gas phase by dust grains. 

The large collecting power of the 10m Keck telescope has made it feasible to study elemental 
abundances in damped Lya systems in much more detail (Wolfe et al. 1994; Lu et al. 1996a). The 
improvements come in three facets. The first is the higher spectral resolution (FWHM< 10 km 
s"-*^), which allows one to see more component structure in the absorption. The second is the much 
improved S/N of the spectra, which makes it considerably easier to study the weak absorption 
lines. Weak absorption lines are crucial for the determination of elemental abundances since they 
do not suffer nearly as much as the strong lines from line saturation. The third has to do with 
the large spectral coverage of the Keck High Resolution Spectrometer (HIRES), which enables 
the study of many elements simultaneously. This is important because the relative abundances of 
elements contain valuable information about the nature of the stellar nucleosynthesis that made 
these elements (cf. Wheeler, Sneden, & Truran 1989). 

For a number of years, we have been collecting high quality spectra of quasars using the Keck 
HIRES in order to carry out a comprehensive analysis of the various absorption phenomena. One 
particular aspect of this program is the study of elemental abundances in damped Lya systems and 
their implications for galactic chemical evolution and for stellar nucleosynthesis. Some preliminary 
results of this nature have been reported in Lu, Sargent, Sz Barlow (1995a). The current paper 
presents the full analyses. The results described here supersede the earlier results reported in Lu 
et al. (1995a). 



-4- 



In §2 we describe the observations, data reduction, and the methods used to estimate 
equivalent widths and column densities. Abundance determinations for the 14 damped Lya 
systems studied here are discussed in §3. These abundance measurements arc combined with 
similar measurements in the literature to study the chemical evolution of damped Lya galaxies 
over the redshift range 0.7 < z < 4.4 (§4). In §5 we discuss other observed properties of damped 
Lya galaxies and implications for their nature. The redshift dependence of the temperature of the 
cosmic microwave background as estimated from the C II* A1335 fine structure line is discussed in 
§6. §7 summarizes the main conclusions. 

Throughout this study, we will assume qo = 0.1 and Hq = 50 km s~^ Mpc~^. We will use z 
for general references of redshift, Zem for the emission redshift of quasars, Zdamp for the redshift 
of damped Lya systems, and Zabs for any other absorption systems or for absorption systems in 
general. 

2. OBSERVATIONS, DATA REDUCTION, AND MEASUREMENTS 

2.1. Observations and Data Reduction 

Table 1 gives the journal of the Keck HIRES observations. A 0.86" slit width was used in all 
cases to yield a resolution of FWHM=6.6 km s~^ with roughly 3 pixels per resolution element. 
The echellc format of HIRES is such that one can only get complete coverage of the free spectral 
range for A < 5100 A in a single setup with the current TeK 2048x2048 CCD. Except in a few 
cases, we used two partially overlapping setups for each wavelength region in order to get complete 
spectral coverage. Data reductions were done using a software written by T. A. Barlow. After the 
echellc orders were optimally extracted and wavelength- and flux-calibrated, they were rcsamplcd 
to 2 pixels per resolution element, scaled to the same flux level, and added together weighted 
according to their S/N. Because some spectral regions have more exposures than others due to the 
overlapping-setup scheme, the S/N of the spectrum from even a single object can show quite large 
variations from one region to another. 

Most metal lines associated with damped Lya systems are blended with Lya absorption lines 
if they occur in the Lya forest. Thus our analyses will primarily rely on measurements made 
redward of the Lya emission line. We estimate the continuum level redward of the Lya emission 
by fitting cubic splines to regions deemed free of absorption lines using the IRAF task continuum. 
In rare cases, we may also use measurements made in the Lya forest if the lines are reasonably 
clean. In such cases, the continuum level is estimated by picking out small "peaky" regions in the 
spectral region of interest that are free of obvious absorption lines and interpolating them with 
low-order cubic splines. The resulting continuum appears to describe the data quite well. 
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2.2. Equivalent Width and Column Density Measurements 

For each damped Lya system studied here, we give a table listing the absorption lines 
detected, the measured equivalent widths, and column densities (see §3). 

To carry out the measurements for a given damped Lya system, we first plot the continuum- 
normalized profiles of the absorption lines in velocity space for an adopted redshift for the system 
(which generally corresponds to the strongest component in the low ion absorption). Equivalent 
width for each line is then estimated by integrating the observed spectrum over the width of the 
absorption line defined by two limiting velocities: V- and v^. The la uncertainty associated with 
the measured equivalent width is also calculated including both the effects of statistical noise and 
of uncertainties in the placement of the continuum level (cf. Savage et al. 1993). Only features 
over the 4(T significance level will be considered a real detection. 

To obtain column densities, we use the apparent optical depth method (cf. Savage & Sembach 
1991). We first convert the normalized absorption line profiles, I{v), into measures of the apparent 
optical depth per unit velocity, Ta{v), through the relation Ta{v) =ln[l//(?;)], which can then be 
converted into the apparent column density per unit velocity, Na{v), in units of atoms cm^^ (km 
s~^)~^ through the relation 

iV„(z;) = 3.767 xlO^^^, (1) 

where A is the rest wavelength in A and / is the oscillator strength. The total column density, 
obtained by integrating Na{v) over the width of the absorption profile, should represent the 
true column density to within the measurement uncertainty if the absorption line contains no 
unresolved, saturated components]^. In the case that the absorption line does contain unresolved, 
saturated components, the integrated Na{v) only gives a lower limit to the true column density. 
For an ion species for which several absorption lines with different fX values are observed, 
comparisons of their integrated Na{v) will generally allow one to infer if unresolved saturated 
components are present. This is because lines with larger /A values will yield lower integrated 
column densities than lines with smaller /A values in such cases. This technique was explored in 
details by Savage & Sembach (1991). Examples of applications may be found in Lu et al. (1995b). 

Empirically, we find that, given the spectral resolution and S/N typical of the data used in 
this work, absorption lines with peak absorption depth less than ~ 80% (i.e., peak optical depth 
To < 1.6 if the lines are resolved) show no evidence of unresolved saturation that would render 
invalid the column densities derived from Na{v) integrations. For example, the Fe II AA2249, 
2260 lines differ in their /A values by a factor of 1.3, and they yield the same A'^(Fe II) for 
several systems in Q 0449—1326, Q 0450—1312, and Q 0454+0356 even though the intrinsically 



•^In principle, one can measure the column density of a saturated absorption line if the absorption line is fully 
resolved, and if the S/N is good enough to clearly indicate that the flux in the center of the absorption line does not 
go to zero. In practice, this rarely happens due to limited S/N in real data. 
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stronger A2260 absorption lines axe 60-75% deep. Similarly, the Si II AA1304, 1526 absorption 
lines, which differ by a factor of 1.5 in their /A values, yield consistent A^(Si II) in several systems 
in Q 1946+7658, Q 2212-1626, and Q 1202-0725 (Lu et al. 1996a) even though the intrinsically 
stronger A1526 absorption lines are 70-85% deep. Also, the S II AA1250, 1253 lines (which differ in 
their /A values by a factor of 2) in the Zdamp = 2.8110 system toward Q 0528—2505 yield consistent 
column densities despite the fact that the intrinsically stronger A1253 absorption is 85% deep. 
The most dramatic example is offered by the Zdamp — 

2.8268 system toward Q 1425+6039, where 
the Fe II AA1608, 1611 absorption lines, which differ in their /A values by a factor of 62, yield 
A^(Fe II) values differing by less than 0.1 dex, and the A1608 absorption line is 80% deep. Similar 
examples may be found for the more ionized species such as Al III, Si IV, and C IV, each having 
doublet absorption lines differing in their /A values by a factor of 2. In fact, we have yet to find an 
example where an absorption line with peak absorption depth < 80% shows convincing evidence 
of being saturated and unresolved. Thus, for ion species for which only a single absorption line is 
observed (i.e., for which we cannot check independently if there is any unresolved saturation in 
the absorption), we will assume that the column density derived from integrating its Na{v) profile 
is unbiased with respect to unresolved saturation if its peak absorption depth is less than 80%. In 
reality, most of the single absorption lines studied here are < 50% deep, corresponding to a peak 
optical depth of 0.7 if the lines are resolved. The only single absorption lines with peak absorption 
depth between 50 to 80% for which we have adopted their A^a(i;)-integrated column density as 
real measurements (rather than lower limits) are the Si II A1808 lines at Zdamp = 2.8110 toward Q 
0528-2505 and at Zdamp = 4.0803 toward Q 2237-0608, the Fc II A1608 line at Zdamp = 2.8268 
toward Q 1425+6039, and Al II A1670 lines at Zdamp = 2.8443 toward Q 1946+7658 and at 
zoamp = 4.0803 toward Q 2237-0608 (see §3 for details). 

We emphasize that, for unsaturated or mildly saturated lines, the apparent optical depth 
method will yield column density measurements in agreement with the profile fitting technique. 
For heavily saturated lines, neither method will yield very reliable column densities in practice. 
For this reason, the analyses in this paper will concentrate on column density measurements made 
from unsaturated lines. Column density limits from saturated lines may be used occasionally in 
the discussion. For lines that arc undetected or detected at significance level less than 4(T, we will 
use the 4cj upper limits of the equivalent width to estimate the upper limits of the column density 
assuming linear part of the curve of growth. 

Atomic data, including rest-frame vacuum wavelength, oscillator strength (/-value), and 
radiation damping constant, are taken from Morton (1991) except for those which have revised 
/-values as given in the compilation of Tripp, Lu, & Savage (1996; see Savage & Sembach 1996 for 
additional information) . 
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3. DISCUSSION OF INDIVIDUAL SYSTEMS 

In this section, we derive column densities and elemental abundances for each individual 
damped Lya system studied here. In general, absorption lines occurring in the Lya forest will not 
be discussed unless the line is relatively clean of contamination from forest absorption. It will be 
assumed that the ionization state of the absorbing gas is dominated by species appropriate for H I 
gas (e.g., O I, N I, C II, Si II, Fc II, etc.) so that abundance corrections for unobserved ionization 
stages are negligible. This is believed to be a good assumption given the large H I column density 
of the systems (see Viegas 1995). 

Two damped Lya systems not covered by the spectra listed in Table 1 will also be discussed. 
These are the z^amp — 

1.1743 system toward Q 0450-1312, and the 

^damp — 0.8598 system toward 
Q 0454+0356. Keck HIRES observations of these two quasars were made by Drs. Christopher 
Churchill and Steven Vogt, who have kindly made their spectra available to us. Details of the 
observations and reductions for these two objects are discussed in Churchill k. Vogt (1997). 

3.1. Q 0000-2620 Zdamp = 3 3901 Damped System 

This damped Lya system has been known ever since the quasar was first discovered by C. 
Hazard (unpublished; but see Webb et al. 1988). It is one of the four damped Lya galaxies at 
z > 3 for which abundance information has been obtained. The HIRES spectrum used here for this 
object was already published by Lu et al. (1996b). We estimate iV(H I)=(2.6 ± 0.5) x 10^^ cm~^ 
from fitting the Lya damping profile (figure 1), which agrees quite well with other determinations 
(Webb et al. 1988; Sargent, Steidel, & Boksenberg 1989; Savaglio, D'Odorico, & Moller 1994). 
Absorption lines detected redward of the Lya emission are listed in Table 2; the absorption profiles 
are shown in figure 2. Several other lines (Si II AA1190, 1193, 1260, 1304, Si III A1206, O I A1302, 
C II A1334, Si IV AA1393,1402) are clearly detected in the Lya forest but are blended with Lya 
forest absorption (see the spectrum in Lu et al. 1996b). These lines will not be discussed here. 

Except for the line tentatively identified as S II A 1253, all the low ion lines shown in figure 2 
are clearly saturated, so only lower limits to their column densities are given in Table 2, along with 
the abundance limits. Similarly, only a lower limit to the C IV column density is given. For Fe II, 
we give both the lower and upper limits to the column density determined from the saturated Fe 
II A1608 absorption and from the absence of the Fe II A1611 absorption. 

In Lu et al. (1996b) we tentatively identified a weak absorption feature in the Lya forest as 
the S II A1253 line associated with this damped Lya system. This weak feature is also shown in 

figure 2, and it shows excellent redshift agreement with the other low ion lines. The corresponding 
S II AA1250,1259 lines are badly blended with Lya forest absorption. If our identification of 
this S II line is correct, then the estimated column density, log A''(S II)=14.70 it 0.03 (Lu et al. 
1996b), would imply an abundance of [S/H]=— 1.98. This is consistent with the limits obtained 
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for the other elements. We note that Pettini et al. (1995b) estimated [Zn/H]< —1.76 (3cr) and 
[Cr/H]= -2.46. 

3.2. Q 0216+0803 Zdamp = 2.2931 System 

This damped Lya system was discovered by Lanzetta et al. (1991) and confirmed by Lu 
& Wolfe (1994), who estimated log A^(H I)= 20.45 ib 0.16 for the system. Our Keck spectrum 
does not cover the wavelength region containing the damped Lya absorption. We will therefore 
adopt the N{}i I) value from Lu Sz Wolfe (1994). The absorption lines detected redward of the 
Lya emission are listed in Table 3, and the absorption line profiles are shown in figure 3. Lines 
detected in the Lya forest, which will not be used in the analysis, include C II A1334, O I A1302, 
Si II A1304, and Si IV AA1393,1402. 

The Zn II A2026 and Cr II A2056 absorption lines are detected at the 2.3cr and 2.6a levels, 
respectively; we will only use their 4a" upper limits in the subsequent analyses. The Fe II A1608 
absorption is probably saturated and we conservatively give a lower limit of logiV(Fe !!)> 14.81. 
The much weaker Fc II A2260 absorption is detected at the 2.6cj level, and the corresponding 4a 
upper limit is log iV(Fe II) < 14.97. We will therefore adopt log iV(Fe II)= 14.89 ± 0.08. 

3.3. Q 0216+0803 Zdamp = 1-7688 System 

This system was identified as a candidate damped Lya absorption by Lanzetta et al. (1991) 
based on a low resolution spectrum. No higher-resolution confirmation of the damped Lya 
absorption is available. However, the fact that we detect intrinsically weak absorption lines like Fe 
II A2260, Mn II A2576, and Ni II A1751 suggests that most likely this is a damped Lya system. 
For example, even for solar metallicity, the observed Fe II column density would imply A^(H 
I)~ 10^^ cm^^. The Lya line equivalent given by Lanzetta et al. (1991) implies log A^(H I)=20. 
We tentatively adopt logA^(H I)= 20 ± 0.18, assuming a 50% uncertainty. The absorption lines 
redward of Lya emission for this system are listed in Table 4, and the profiles are shown in figure 
4. Note that none of the Zn II and Cr II absorption lines is detected at over 4(7 significance level 
so only upper limits to their column densities are listed in Table 4. Metal lines detected in the 
Lya forest include C IV AA1548,1550, Fe II A1608, and Al II A1670. 

3.4. Q 0449-1326 Zdamp = 1 2667 System 

The redshift of this system is too low to see the corresponding Lya absorption line {Xobs = 2755 
A) from the ground. The damped Lya nature of the system is inferred from the detection of the 
intrinsically weak Mn II triplet lines and the Fe II AA2249,2260 lines (see Table 5), which even for 
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solar abundances, implies a A/"(H 1)^ 4 x 10^^ cm ^. The metal lines detected in this system are 
listed in Table 5. Figure 5 shows the absorption line profiles. 

The derivation of the Mn II column density deserves some discussion. The estimated values 
of A'^(Mn II) obtained by integrating Na{v) increases systematically with decreasing /A value of 
the lines (Tabic 5). This seems to indicate that all three Mn II lines are moderately saturated. 
However, we argue that this is most likely not the case. The Mn II lines are weaker than or 
comparable to the Fe II A2260 absorption in the system. The good agreement in the estimated 
N{Fe II) from the Fe II A2260 absorption and the weaker Fe II A2249 absorption indicates that 
the Fe II A2260 absorption is not saturated. Thus it is difficult to understand why the Mn II lines 
should be saturated. The S/N near Mn II A2576 is twice as good as the region near Mn II A2594, 
and three times as good as the region near Mn II A2606. The values of A^(Mn II) from the A2576 
and A2594 lines agree well within la. It is only the iV(Mn II) from the A2606 line that differs 
significantly. A careful look at the actual apparent column density profiles (figure 6) of the three 
Mn II lines confirms that it is the Mn II A2606 absorption that shows significant deviation from 
the other two lines. It is likely that the Mn II A2606 absorption is contaminated by an unidentified 
weak absorption line from another redshift system. We will therefore adopt the Mn II column 
density determined from the A2576 and A2594 lines. 

Since we do not have the A'^(H I) information, it is not possible to derive absolute abundances 
for the system. However, it may be interesting to get a crude idea of what the abundances might 
be for reasonable guess of the A^(H I). For A^(H I)=2 x lO^'^ cm~^, the column densities given in 
Table 5 implies: [Fe/H]=-0.68, [Mn/H]=-0.96, and [Mg/H]> -1.78. As we will see later in §4, 
the overabundance of Fe relative to Mn in this system, [Fe/Mn]=+0.28, is quite representative of 
damped Lya systems in general and has important implications for the nucleosynthetic history of 
the absorbing galaxies. 

3.5. Q 0450-1312 Zdamp = 1 1743 System 

The HIRES spectrum of this object came from Churchill & Vogt (1997). Figure 7 shows the 
metal line profiles, where the intrinsically weak Fe II AA2249, 2260, Mn II AA2576, 2594, 2606, and 
Cr II A2056 lines are clearly detected. The equivalent widths and column densities of the lines are 
given in Table 6. The A^(H I) of the system is unknown so it is not possible to derive absolute 
metal abundances. However, the Fe II column density implies N(H. I)=4 x 10^^ cm~^ even for 
solar metallicity. Hence we believe that this is a damped Lja system. 

3.6. Q 0454+0356 Zdamp = 0.8598 System 

Abundances of Fe, Zn, and Cr for this damped Lya system have been determined by Steidel 
et al. (1995) using a medium-resolution spectrum; they found [Fe/H]= —1.04, [Zn/H]= —1.08, and 
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[Cr/H]= —1.07 (note that these values have been corrected for the /-values of Tripp et al. 1996). 
Steidel et al. also estimated a A'"(H I)=(5.7 ± 0.3) x 10^'^ cm~^ from the damped Lya absorption 
in a HST FOS spectrum of the quasar. The Keck HIRES spectrum used in this study came from 
Churchill &; Vogt (1997). Figure 8 shows the metal lines, while Table 7 gives the measurements. 
The stronger Zn II A2026 line is not covered by the HIRES spectrum. In any case, the S/N of the 
spectrum near the Zn II and Cr II lines is very low, and none of the Zn II or Cr II lines is detected 
at > 4a" significance level. The upper limits on [Zn/H] and [Cr/H] (Table 7) are consistent with 
the Stcidcl ct al. (1995) measurements. The [Fe/H] estimated from the HIRES spectrum also 
agrees with the Steidel et al. estimate. The only additional information obtained for this system 
is the Mn abundance. 



3.7. Q 0528-2505 Zdamp = 2.8110 System 

We detect over 30 metal lines from 15 different species in this system, partially because 
the damped Lya system has a redshift slightly higher than the emission redshift of the quasar 
{zem = 2.779), so more metal lines than usual are shifted out of the Lya forest. Even though our 
spectrum covers most of the damped Lya absorption line, the coverage to the blue side of the 
damped Lya absorption is too short to reliably assess the continuum level before the damped Lya 
absorption. We therefore adopt log A^(H I)= 21.2 it 0.1 from Morton ct al. (1980) and Foltz, 
Chaffee, & Black (1988). The metal lines are presented in Table 8 and figure 9. We also derive 
upper limits to the abundance of several other elements, including P, CI, Ti, Co, Cu, Ga, and Ge, 
in order to compare with the results of Songaila & Cowie (1996). Note that our abundance upper 
limits are 4cr values. The corresponding 2a limits should be 0.3 dex lower than those given in the 
table. 

This system is unique in several ways. First, the redshift of the damped Lya absorption is 
higher than that of the quasar, suggesting that the absorbing galaxy may be physically associated 
with the quasar in a group or cluster. The presumed close proximity of the background quasar, 
which is an intense source of ionizing photons, may have significant effect on the physical state 
of the absorbing gas in the damped Lya galaxy. Secondly, this is the only damped Lya system 
for which molecular hydrogen absorption has been detected (Foltz et al. 1988; Songaila &; Cowie 
1996). Thirdly, the metal absorption lines are unusually wide and complex, spreading over 600 
km s^^ in the strongest low ion lines and 500 km s~^ in the high ion lines. It appears that the 
absorption consists of two groups of components at above and below 140 km s^"*^, respectively 
(figure 9). The complexity of the metal line absorption may be caused by two galaxies in 
superposition. 

The Si II A1808 absorption is relatively strong (about 70% deep) and might be somewhat 
saturated. However, the comparably strong S II A1253 absorption does not appear to show any 
unresolved saturation based on a comparison of the N{S II) column densities derived from this 
line and from the weaker A1250 line (Table 8). We will therefore assume that the N{Si II) derived 
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from the A 1808 absorption is correct. 

Of particular interest is the C II* A1335 absorption. The ratio of iV(C ll*)/N{C II) is a density 
indicator (§5). The ratio can also be used to estimate the temperature of the cosmic microwave 
background radiation at the redshift of the absorption system (§6). The C II* absorption in this 
system is clearly detected but is partially blended with the C II A1334 absorption line (figure 9). 
However, we note that the part of the C II* absorption at v > 140 km s~^ is free of any obvious 
blending. Hence it is possible to estimate N{C II*) for this part of the absorption (see last section 
of Table 8). 

Another ion of special interest is N I since the abundance ratio of N/0 can be used to infer the 
origin of the heavy elements (§4). We have clearly detected N I absorption from the triplet near 
1200 A, but they are all blended with each other (figure 9). However, the part of the absorption 
Sit V > 150 km s-i for the weakest N I line A1200.710 appears to be free of blending. Thus we 
have attempted to derive A^(N I) for this part of the absorption, and also for a few other elements 
(see last section of Table 8) in order to examine the relative abundances. A direct integration of 
the apparent column density profile of the N I A1200.710 line between v = [+140, +300] km s~^ 
yields log N{N I)= 14.86 ± 0.02. It is possible that the N I A1200.710 absorption in this velocity 
interval is somewhat saturated. In order to be conservative, we use the Fe II lines to estimate 
the maximum possible degree of saturation in the N I absorption. The Fe II AA1608,2260 lines 
yield log A^(Fe II)= 14.65 ± 0.02 and log A^(Fe II)< 15.12 (4(7) respectively, over the interval 
V = [+140, +300] km s~^ (Table 8). The Fe II A1608 absorption might be saturated in this velocity 
interval, but the column density from the A2260 line indicates that the value of log A^(Fe II)=14.65 
may underestimate the true N{Fe II) by at most 0.47 dex. Since the N I A1200.710 absorption 
has similar optical depth as the Fe II A1608 absorption over the velocity interval [+140, +300] km 
s~^, it may be reasonable to assume that the true log A^(N I) over this velocity interval should be 
no more than 14.86+0.47=15.33. Thus we find conservatively log 14.65 < A^(N I)< 15.33, and 
[N/0]< +0.65 or -1.21 <[N/Si]< -0.53 for the absorption complex at v > 140 km s~^. 

3.8. Q 0528-2505 Zdamp = 2.1410 System 

Morton et al. (1980) estimated A'^(H I)= 5 x lO^'' cm~^ for this system, which we will 
adopt and assume a 20% uncertainty. Figure 10 shows the metal absorption lines detected in 
this system. The measurements are given in Table 9. The Zn II A2026 absorption is probably 
detected but unfortunately blended with the Al II A1670 line at z^amp = 2.8110, and we only 
have an upper limit to N(Zn II) from the weaker Zn II A2062 line. The Fe II AA2249,2260 lines, 
which normally provide the best measure of A^(Fe II) owing to their weak but detectable strength, 
are unfortunately blended with the Al III A1854 and A1862 lines in the Zdamp = 2.8110 damped 
system, so no measurements are possible. We can only derive a lower limit to A''(Fc II) from the 
moderately saturated Fe II A1608 line. An upper limit of log N{Fe II)< 15.21 can be obtained 
from the absence of Fe II A1611 absorption. We thus adopt log N(Fe II)=14.94 ± 0.26 in the 
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analysis. 

3.9. Q 1425+6039 Zdamp = 2.8268 System 

This bright high redshift quasar {zem = 3.2, V=16.5) was discovered by Stepanian et al. 
(1991). Our Keck spectrum for this object is of exceptionally high quality, with typical S/N 
reaching 120 per resolution element redward of Lya emission. The damped Lya absorption line 
in the Keck spectrum is shown in figure 11. The metal absorption lines are shown in figure 13. 
A number of lines are detected in the Lja forest, including Si II AA1190, 1193, 1260, 1304, Si III 
A1206, O I A1302, and several Fe II lines. These will not be used in the analysis. 

The metal absorption lines in this damped system contain two groups of components, one 
between [-100, +100] km s~^, and another one between [+200, +300] km (figure 13). The 
low ion absorption lines clearly indicate that the bulk of neutral gas occurs at zero velocity 
{z = 2.82680). However, damping profiles fitted to the Lya and Ly/3 absorption lines, with the 
redshift fixed at 2.82680, yield satisfactory fits to the blue sides of the profiles for A^(H I)=2 x lO^'^ 
cm~^ and 6 = 10 km s~^, but leave a big chunk of absorption unaccounted for at the red side of 
the line profiles (figure 11). The "extra" absorption is clearly due to the group of components 
occurring between +200 and +300 km s~^. Figure 12 shows the same as figure 11, but with 
another component added at Zabs = 2.83058, corresponds to v = 300 km s~^ in figure 13. The 
extra component has A^(H I)=1.0 x 10^^ cm~^ and 6 = 10 km s~^, which fits the Lya profile well, 
but not that of Ly/?. This is because the red side of the Lya profile is dominated by the damping 
wing of the component at v = 300 km s~^ (i.e., independent of the velocity dispersion of the gas), 
while the Ly/3 line, saturated but not damped, depends on the relative distribution of the gas in 
the group of components between [+200, +300] km s~^ and on the velocity dispersion of each 
component. Figure 12 suggests that the total A'^(H I) of the group of components between [+200, 
+300] km s~^ is roughly 10"*^^ cm^^. Experiments indicate that the total A^(H I) of this group is 
definitely no higher than 2 x 10^*^ cm^^. Because of this relatively low A^(H I), one might worry 
about ionization corrections. For this reason, we will confine the analysis of this system to the 
components at v < 120 km s"^, and adopt a log A^(H I)= 20.3 ± 0.03. Table 10 lists the metal 
lines and the measurements. 

The N I triplet near 1200 A is clearly detected (see figure 13), although the N I A1200.710 line 
(not shown) is contaminated by Lya forest absorption. We are only able to derive a lower limit 
to N^N I). The O I A1302 absorption occurs in the Lya forest and is strongly saturated. Thus 
even deriving a limit to the N/0 ratio is not possible. Our spectrum docs not extend red enough 
to cover the weak Si II A1808 absorption line. The Si II A1526 absorption is saturated, and the 
component near —90 km s~^ is blended with the C IV A1548 absorption at Zabs = 2.7727. We are 
only able to derive a lower limit to A'^(Si II) from the Si II A1526 absorption at v > —50 km s"^. 
The Fe II A1608 absorption is relatively strong and might be somewhat saturated. We estimate 
log iV(Fe II)= 14.44 from integrating its Na{v) profile. The much weaker Fe II A1611 absorption is 
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detected with Wr = 0.007 ±0.002 A (3.5cr detection), with a column density of log N{Fe II)=14.48. 
This suggests that the Fe II A1608 absorption is unsaturated or nearly unsaturated. However, in 
order to be self-consistent, we adopt log A^(Fe II) < 14.52 {4a) for Fe II A1611, and adopt log A^(Fe 
II)= 14.48 lb 0.04 as the final column density. The S II lines arc apparently blended with strong 
Lya absorption lines so no measurements are possible. Most (perhaps all) of the absorption in the 
C II* A1335 panel is from the C II A1334 line. The contamination made it impossible to estimate 
the N{C II*) or its upper limit. 

3.10. Q 1946+7658 Zdamp = 2.8443 System 

This damped Lya system and another one at Zdamp = 1.7382 in the same spectrum (see 
§3.11) have been studied in some detail by Lu et al. (1995b; see also Fan & Tytler 1994). The 
Keck spectrum obtained here is of higher quality than those used in previous studies and extends 
further into the blue. Thus it is worthwhile to repeat the analysis, particularly in light of the N I 
lines which are covered by the Keck spectrum. The A^(II I) of the system is estimated to be log 
A?^(H I)= 20.3 lb 0.1 by Lu et al., which agrees well with the new determination based on the Keck 
spectrum: log N{E. I)=20.27 ib 0.06 (see figure 14). 

The metal lines in this systems are shown in figure 15. The Si II A1190,1193 and Si III A1206 
lines are detected in the Lya forest. The comparison between the C II A1334 absorption and 
the O I A1302 absorption is very informative. While the O I absorption has only one obvious 
component centered at zero velocity, the C II absorption shows at least two extra components at 
—70 and +80 km s~^. This difference is almost certainly not due to the relative strength of the O 
I and C II absorption since the much weaker Al II A1670 absorption (and possibly the Si II A1526 
absorption) also shows the two extra components as seen in C II. The difference is most likely due 
to ionization effects. The O I ion may be considered the best tracer of neutral gas (other than H I 
itself) since it has an ionization potential (13.61 eV) nearly identical to that of H I. On the other 
hand, ions like C II, Al II, Si II, etc. can exist in partially ionized regions because of their higher 
ionization potentials. Thus essentially all of the neutral hydrogen in the system is associated with 
the component at zero velocity. We will hence restrict the column density measurements of the 
low ion lines to the component at zero velocity (which is reflected in the v- and values used 
to estimate the equivalent width and column density). This distinction makes an appreciable 
difference only for Al II. The results are given in Table 11. 

The C II* A1335 absorption is clearly absent, and the upper limit obtained here is twice as 
good as that given by Lu et al. (1995b). By chance, all three N I lines occur in clean regions of the 
Lya forest (figure 15; for space considerations the weakest N I line is not shown); however, none 
is detected at > 4c7 significance level. The upper limit on A^(N I) and the lower limit on A'^(0 I) 
from the saturated O I A1302 absorption yields a firm upper limit of [N/0]< —0.76 dex. We are 
also able to put a not- very-informative upper limit on [S/H] from the absence of the S II A1250 
absorption, which happens to occur in a clean region of the forest. The two stronger S II lines are 
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blended with Lja absorption. The absorption in the N V A 1242 panel is most Hkely due to Lya 
forest hnes. 



3.11. Q 1946+7658 Zdamp = 1-7382 System 

The Lya Une at this redshift occurs at a wavelength blueward of the Lyman limit absorption 
from the z^amp = 2.8443 damped system in the same quasar spectrum and is thus not observable. 
The damped Lya nature of the system is inferred from the detection of the intrinsically weak Si 
II A1808, Fe II AA2249,2260 and the Cr II lines (figure 16), which implies N{R I)= 1.6 x 10^^ 
cm~^ even for solar abundances. The metal lines detected in the system are shown in figure 16 
and listed in Tabic 12. The Cr II A2066 absorption is blended with the C IV A 1548 absorption at 
Zabs = 2.6541. The C IV absorption lines in this damped system occur in the Lya forest and is 
probably somewhat contaminated. 

The column densities for Al III obtained from the Al III AA1854,1862 absorption seem to 
indicate unresolved saturation in the Al III lines as the difference between the two estimates, 
0.12 dex, is significantly larger than the la errors. This is surprising given the relative weakness 
of the absorption lines. However, most of the difference occurs near v = 20 km s~^, where the 
intrinsically weaker Al III A1862 absorption is even deeper than the intrinsically stronger A1854 
absorption. We believe the Al III A1862 absorption is probably contaminated by an unidentified 
metal absorption line. We will assume this is the case and adopt the N{Al III) from the A1854 
line. 

The lack of N(H I) information prevents an estimate of the absolute abundances in this 
system. For reference, if log A^(H I)=20.3, then the following abundances would be obtained: 
[Si/H]= -1.09, [Fe/H]= -1.35, [Zn/H]=< -1.12, [Cr/H]= -1.20, and [Mn/H]= -1.71 (adopting 
log iV(Mn II)=12.12 from Lu et al. 1995b). 

3.12. Q 2212-1626 Zdamp = 3.6617 System 

This is the third highest redshift damped Lya system for which detailed analysis has been 
carried out. Only the Zdamp — 

4.0803 damped system toward Q 2237-0608 (§3.14) and the 
^damp = 4.3829 damped system toward BR 1202—07 (Lu et al. 1996a) axe at higher redshifts. 
From our Keck spectrum (see figure 17), we estimate log A^(H I)=20.20 ± 0.08. The metal lines 
are shown in figure 18 and listed in Table 13. A number of other metal lines (e.g., Si III A1206, 
Si II AA1190, 1193, 1260) appear to be present in the Lya forest but are badly blended with forest 
absorption. 

One peculiarity about this system is that the O I A1302 absorption clearly does not reach 
zero flux in its line center, which is unusual for this normally-strong absorption line. The O I 
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absorption occurs on top of the strong Lya emission line of this quasar. Since heavily saturated 
Lya absorption lines to the blue of the O I absorption clearly goes to zero flux in their line 
centers, the residual flux in the O I absorption line center is not an artifact of the data reduction 
procedure. Indeed, visual inspections of the raw images clearly confirm that the flux in the O I 
absorption line center is not zero. We consider two possible causes of this residual flux: resolution 
smearing and ionization effects. 

If the O I absorption consists of many extremely narrow components that are not resolvable 
at our resolution (FWHM=6.6 km s^^), smearing of the components due to the inadequate 
resolution could lift the flux in the absorption line center above the zero level even if the absorption 
components are intrinsically saturated. This requires the components to have b < A km s^^ or 
temperature T < 1.5 x 10^ K. Experiments show that, indeed, one can explain the residual flux 
in the O I absorption line center if it contains several appropriately spaced components narrower 
than 6 = 4 km s~^. 

Another possible explanation of the residual flux in the O I absorption is that it is real, 
namely, the O I absorption line is unsaturated. This is unusual but not impossible. If this is 
the case, then the A^(0 I) from integrating the Na{v) profile (see table 13) should be an actual 
measurement rather than lower limit. This in turn implies [0/H]=— 2.37 for this system. This O 
abundance is nearly 0.5 dex below the Si abundance ([Si/H]=— 1.90) in the same system, which 
is inconsistent with nucleosynthesis considerations (see §4.3). One can reconcile this difference if 
the Si abundance in this system has been overestimated because some of the Si II ions come from 
partially ionized gas rather than from neutral gas^ If this is the case, it would mean substantial 
ionization corrections for this system. However, we consider this explanation unlikely for the 
following reasons. Firstly, when there are significant ionization corrections, usually there are 
notable differences in the absorption profiles between O I and other ions that can exist in H II 
regions (e.g., C II, Si II, etc.; see the Zdamp = 2.8443 system toward Q 1946+7658 in §3.10). In 
contrast, the O I absorption in the Q 2212—2616 system has identical velocity structure to the C II 
and Si II absorption in the system. Secondly, the upper limit on the temperature of the microwave 
background radiation at the redshift of this absorption system as derived from the C II*/CII 
ratio (§6) is extremely close to the Big Bang prediction, indicating that collisional excitation of 
the C II atoms by free electrons is negligible in the absorbing gas. This is in contradiction with 
the large electron density expected in this system if the absorbing gas is substantially ionized. 
These considerations suggest that the resolution-smearing interpretation given in the preceding 
paragraph is likely the correct explanation. In fact, the implied low electron density in this system 
from the C II*/CII ratio (§6) suggests that this system may be less ionized than typical damped 
Lya systems, which in turn may imply a lower temperature and hence narrower line widths (see 
the preceding paragraph). 



^Recall that it takes 16.3 eV to turn Si II into Si III, so Si II ions can exist in H II regions. In contrast, the 
ionization potential ol O I is 13.61 eV, making it the best tracer of neutral gas other than H I itself. 
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3.13. Q 2231-0015 Zdamp = 2.0662 System 

Our Keck spectrum does not cover the damped Lya line. Lu &: Wolfe (1994) estimated log 
A'^(H I)=20.56 ± 0.10 for this system, which we adopt. The metal lines are shown in figure 19 and 
listed in Table 14. Metal lines detected in the Lya forest include Si II A1526, C II A1334, Si IV 
AA1393,1402, and C IV AA1548,1550, but they all appear to be contaminated by forest absorption. 
The Fc II A1608 absorption yields a log A^(Fc 11)= 14.72 it 0.01 from integrating the apparent 
optical depth profile. Given its strength, it is probably saturated. Thus we take the above value 
as lower limit. The absence of the Fe II A1611 absorption yields an Aa upper limit of log A'^(Fe 
II)< 15.08. Hence we conservatively adopt log iV(Fe II)= 14.90 ± 0.18. The Cr II absorption is 
below the 4(7 significance level. 

3.14. Q 2237-0608 Zdamp = 4.0803 System 

This is the second highest redshift damped Lya galaxy studied in detail. The damped Lya 

absorption line in the Keck spectrum is shown in figure 20. We estimate a iV(H I)=3 x 10^° cm~^ 
with about 30% uncertainty. This A^(H I) estimate is more uncertain than the estimates for the 
other damped systems based on Keck spectra because both wings of the damped Lya line are 
affected by strong, dense Lya forest absorption due to its very high redshift. The iV(H I) may 
have been slightly underestimated judging from the base of the damped Lya absorption, although 
increasing N{H. I) will overfit the wings farther away from the line center. 

The metal lines in this system are shown in figure 21; the measurements are given in Table 
15. The C II A1334 and C II* A1335 lines (and to some degree the Ni II A1370 line) happen to 

occur in the red wing of the strong Lya emission line so the S/N are exceptionally good. The Si 
II A1526 and Al II A1670 absorption lines are about 80% and 70% deep, respectively. Based on 
the discussion in §2.2, they arc likely to be unsaturated or nearly unsaturated. We will therefore 
adopt their column densities as measurements rather than lower limits. 

Many metal lines in this system are expected in the Lya forest, but trying to identify them is 
meaningless given the strong possibility of contamination from the extremely dense Lya forest. 

4. CHEMICAL EVOLUTION OF DAMPED LYMAN-ALPHA GALAXIES 

In this section we combine the abundance information obtained for the sample of damped 
Lya galaxies discussed in §3 with those from previous studies in order to carry out a statistical 
analysis of the chemical evolution of damped Lya galaxies. In particular, we will examine the 
metallicity distribution of damped Lya galaxies and compare that with those for the various 
components (disk, halo) in the Milky Way galaxy and with external galaxies. We will examine the 
age-metallicity relation of the damped Lya galaxies and compare that with the similar relation 
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determined for the Milky Way. We will also examine the abundance ratios of various elements 

with different nuclcosynthetic origin in order to gain more insight into the chemical enrichment 
processes in these damped Lya galaxies. The effects of dust on the derived abundances through 
selective depletion will also be discussed. 

4.1. Sample Construction 

4.1.1. Damped Lya Abundances 

Elemental abundances have been measured for a number of damped Lya galaxies in the 
literature with varying degrees of accuracy. In constructing the sample of damped Lya galaxies 
for further statistical analysis, we will only consider those abundance measurements that we 
believe to be largely free of saturation effects. This means that, in general, measurements (but 
not lower limits) from heavily saturated lines will be discarded. Previous abundance estimates 
or upper limits that were based on lines with less than 4cr significance level will be converted 
where possible into 4(7 upper limits to conform to the standard adopted in §3, or else they will be 
discarded. Where appropriate, we will also correct the published abundance measurements for the 
set of new oscillator strengths compiled by Tripp ct al. (1996) so that all the measurements will 
be on the same footing. Following tradition, we will use Fe as the metallicity indicator. Hence, in 
general, only damped Lya galaxies whose Fe abundance is available will be included in the sample. 
However, as we will see in §4.3, the Cr abundance tracks Fe abundance very well in damped Lya 
galaxies, i.e., the Cr/Fe abundance ratios in damped Lya galaxies (when both are measured) are 
consistent with their solar ratio. In order to increase the sample size, we will sometimes substitute 
[Cr/H] for [Fe/H] when the Fe abundance is not available (this happens in only 4 of the 23 systems 
listed in Table 16). Hence damped Lya galaxies without [Fe/H] measurement, but with [Cr/H] 
measurement, will be included in the sample. The final sample of damped Lya galaxies, including 
those discussed in §3 and those from previous studies that satisfy the above selection criteria, is 
given in Table 16. 

4.1.2. Galactic Disk and Halo Star Abundances 

The intrinsic chemical compositions of Galactic disk and halo stars of different age reflect 
the conditions of the local interstellar medium in the Galaxy when the stars were formed, which 
depend on the cumulative effects of the rate of star formation, the stellar initial mass function, 
stellar yields, gas infall/outflow, etc. prior to that epoch. The chemical compositions of Galactic 
stars not only allow us to study the past history of chemical evolution of the Milky Way, but 
also provide a fundamental reference of what might be expected of the chemical compositions of 
interstellar gas at different stages of a galaxy's enrichment history. In §4.3 we will compare the 
relative abundances of various elements found in damped Lya galaxies to those found in Galactic 
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stars in order to gain some insight into the chemical enrichment processes occurring in damped 
Lya galaxies. In constructing the sample of stellar abundances, the following sources are used: 
Gratton & Sneden 1988, 1991 (Fe, Si, Cr, Ni), Edvardsson et al. 1993 (Fe, Si, Ni), Clegg, Lambert, 
& Tomkin 1981 (Fe, S), Francois 1987, 1988 (Fe, S). Gratton 1989 (Fe, Mn), Beynon 1978 (Fe, 
Mn), Sneden k Crocker 1988 (Fe, Zn), Sneden, Gratton, Crocker 1991 (Fe, Zn), and Magain 1989 
(Fe, Cr). 

We note that stellar abundances of sulphur are rarely studied mainly because the most useful 
S I lines occur in the near-infrared (~8700 A) and are weak and blended. Of the three references 
for sulphur that we can find, Clegg et al. (1981) studied 20 stars mostly with [Fe/H]> —0.5, while 
Francois (1987, 1988) contain 26 stars mostly with —1.5 <[Fe/H]< —0.5. Although the Francois 
sample stars are more useful for comparison with the damped Lya galaxies from the metallicity 
point of view, Lambert (1989) argued that the sulphur abundances of Francois for halo stars, 
([S/Fe]) = 0.6 dex, may be too high by 0.2 dex owing to the S I /-values used. Fortunately, this 
uncertainty does not affect any of the conclusions reached later. 

4.1-3. Galactic ISM Cloud Abundances 

It is well known that the relative abundances of elements in Galactic ISM clouds differ 
significantly from the solar composition and vary widely from sight line to sight line, i.e., from 
cloud to cloud. These deviations from solar relative abundances stem from the fact that some 
elements (refractory elements) are preferentially incorporated into solid forms or dust grains, while 
other elements are much less affected. Assuming that the intrinsic abundances of the elements in 
the ISM clouds are the same as the solar abundances, the measured gas-phase abundances of the 
elements in ISM clouds then allow one to infer the fractions of different elements that are locked 
up in grains, i.e., the amount of depletion. Thus the abundance ratios of elements measured in 
diffuse ISM clouds provide a fundamental reference to gauge the effects of dust depletion. With 
the understanding that the properties of dust grains in damped Lya galaxies may be significantly 
different from those of Galactic dust^, we will compare in §4.3 the abundance ratios found in 
damped Lya galaxies with those found in the ISM clouds in order to get a handle on the possible 
presence or absence of dust in damped Lya galaxies. 

Abundances of Si, S, Cr, Mn, Fe, Ni, and Zn measured for ISM clouds based on recent 
high quality HST GHRS data are compiled in Table 17. Only abundance measurements with 
relatively small uncertainties (typically less than 0.1 dex) that are judged to be free of significant 
saturation effects will be used. As with the damped Lya abundances, we have corrected the GHRS 
measurements for the set of new /-values compiled by Tripp et al. (1996) where appropriate. 



^Indeed, the extinction curves found for the LMC and the SMC are known to be different from each other and 
from that found in the Milky Way, indicating that the dust grains in all three galaxies are different. Similarly, the 
reddening law has been found to depend on radius in M31 (lye & Richter 1985). 
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4.2. Age-metallicity Relation 

Figure 22 shows the distribution of [Fe/H] as a function of rcdshift for the damped Lya 
galaxies in our sample (solid circles). This is equivalent to the agc-mctallicity relation commonly 
known in Galactic chemical evolution studies. We have used [Fe/H] as the metallicity indicator, as 

is conventional. In a few cases, [Cr/H] is used in place of [Fe/H] when the latter is not available 
(see §4.3 for justification). Also displayed in figure 22 is the age-metallicity relation for the sample 
of Galactic disk stars ("+" symbols) studied by Edvardsson et al. (1993). The conversion between 
age (starting from the Big Bang) and redshift assumes qo = 0.1 and Hq = 50 km s^^ Mpc~^, which 
gives a current age of the universe of 16.6 Gyrs, consistent with the stellar age determinations. 

As is clear from figure 22, the damped Lya galaxies have Fe-metallicities in the range of 1/10 
to 1/300 solar, or [Fe/H] between —1.0 and —2.5. Hence they represent a population of young 
galaxies that are still in the early stages of their chemical evolution. In fact, as we argue in §4.3, 
the star formation process in most of these galaxies probably started less than 1 billion years 
before the epoch corresponding to their redshift. The iV(H I)-weighted mean metallicity of the 
damped Lya galaxies between 2 < 2; < 3 is [Fe/H]=— 1.56, which may be regarded as the cosmic 
metallicity at {z) = 2.5. 

It is important to note that the true Fe abundances of damped Lya systems could be 
systematically higher than those indicated in figure 22 if significant depletion of Fe onto dust 
grains has occurred. However, as we argue in §4.3, the amount of Fe depletion in the damped Lya 
galaxies in our sample is probably negligible, and should be on average no more than 0.4 dex. 
Hence modest depletion of Fe will not change the qualitative picture presented in this work. 

4-2.1. The Large Spread in [Fe/H] 

As already noted by Pettini et al. (1994), the chemical enrichment processes in damped 
Lya galaxies appear quite inhomogeneous. For example, the [Zn/H] distribution in the sample of 
damped systems studied by Pettini et al. spans at least a factor of 20 or 1.3 dex. The [Fe/H] in 
our sample at z ~ 2 — 3 spans roughly the same range. A metallicity gradient similar to those 
found in the Milky Way and in local spiral disks can only account for a factor of ~2 spread in 
the observed metallicities of damped Lya galaxies (Pettini et al. 1994). If all galaxies in the 
sample formed at roughly the same time, then the large spread in [Fe/H] may indicate widely 
different star formation histories in these galaxies. However, we consider it more likely that the 
large metallicity spread is due to the different formation epoch of the sample galaxies and/or 
that the sample contains a heterogeneous mix of galaxy types. Steidel (1995) found that galaxies 
responsible for the Mgll absorption systems at 2; < 1.6 in quasar spectra come with a variety 
of morphological types, colors, and luminosities. This may explain the large metallicity spread 
in damped Lya galaxies if they sample the same parent population of galaxies. We point out 
that varying amounts of Fe depletion due to dust grains are unlikely to account for the bulk of 
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metallicity spread since [Zn/H] shows similar spread (Pettini et al. 1994), and Zn is not expected 
to be depleted heavily by dust grains (see discussion in §4.3). 

4-2.2. Evidence for z >3 as the Epoch of Galaxy Formation 

There is evidence in figure 22 that the mean metallicity of the sample increases with time. 
This is perhaps not unexpected, but nonetheless reaffirms the notion that heavy metals gradually 
build up with time in galaxies through successive generations of star formation. There was a hint 
of such a trend in earlier studies (Pettini et al. 1995b), but the trend is made much clearer here 
by the addition of several new measurements at z > 3 and by the fact that most of our data 
points are actual measurements rather than upper limits (as is the case in Pettini et al. 1995). In 
particular, all four of the damped Lya galaxies at z > 3 have [Fe/H]< —1.7 (or, [Fe/H]< —2.0 if 
we use the 2a upper limit for the Q2212— 1626 system). In comparison, many of the galaxies at 
2 < z < 3 have reached ten times higher metallicity. Thus the time around z ~ 3 could be the 
epoch of galaxy formation in the sense that it signifies the onset of significant star formation in 
galaxies. Several other lines of evidence appear to point to the same conclusion, as we discuss 
below. 

The first such evidence comes from the space density distribution of quasars. It is found that 
the space density of quasars peaks at 2 < z < 3 and starts to decline at z > 3 (Schmidt, Schneider, 
& Gunn 1995; Warren, Hewett, & Osmcr 1994). By z ~ 4.3, the space density of quasars has 
dropped by a factor of 7 compared to their peak density at z = 2 — 3 (Kennefick, Djorgovski, & de 
Carvalho 1995). If this drop in quasar space density at z > 3 signifies the early phase of structure 
formation from the highest density peaks of the primordial fluctuation, it is perhaps "natural" to 
expect galaxy formation and subsequent star formation to occur at slightly lower redshifts. The 
second supporting evidence comes from studies of the neutral gas content in damped Lya systems. 
Previous studies found that 0,(iampj the cosmological density contained in the H I gas in damped 
Lya galaxies, decreases from z ~ 3 to z = 0, which was interpreted to indicate the transformation 
of gas into stars (Lanzetta et al. 1995; Wolfe et al. 1996; but also see Pei & Fall 1995). Recent 
extension of such studies to even higher redshift reveals that ^damp declines rapidly at z > 3 
(Storrie-Lombardi Sz Wolfe 1996). This may indicate that damped Lya galaxies are still forming 
at z > 3. The third supporting evidence comes from the morphology of high redshift galaxies. 
HST images of z > 3 galaxies identified from their redshifted Lyman continuum break (Steidel, 
Pettini, & Hamilton 1996a; Steidel et al. 1996b; Giavalisco, Steidel, & Macchctto 1996) reveal 
morphological structures that could be interpreted to indicate star formation in the spheroidal 
component of massive galaxies. 

Even if none of these arguments is convincing by itself, taken together, they do consistently 
point to a picture that galaxies started forming at z around or slightly higher than 3. 

Finally, it may be significant that the lowest Fe-metallicity we find for the sample of damped 
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Lya galaxies is about —2.4, which, coincidentally, also represents the lowest metallicity found 

for Galactic globular clusters. Recent observations of extremely low metallicity halo stars 
([Fe/H]< —2.5) reveal that the Galaxy underwent a distinct phase of nucleosynthesis before 
reaching [Fe/H]~ —2.4, characterized by some unusual elemental abundance ratios and large 
variations (Mc Williams et al. 1995; Ryan, Norris, Sz Beers 1996). Possibly, this marks a time 
before which the chemical composition of the gas is only affected by a few Type II supernovae 
(Audouze & Silk 1995), namely, the epoch of formation of the first generation stars. More studies 
of z > 3 damped Lya galaxies may eventually reveal such systems (i.e., [Fe/H]< —2.5), which 
would provide unambiguous evidence for the initial formation of galaxies. 

4-2.3. Comparison with the Milky Way Galaxy and Other Nearby Galaxies 

We note from figure 22 that the degree of chemical enrichment in the sample of damped 
Lya galaxies is, on average, considerably lower than the Milky Way disk at any given time in 
the past. This may bear significantly on the nature of the damped Lya galaxies. It was initially 
suggested (cf. Wolfe 1988) that the damped Lya absorbers may trace disks or proto-disks of 
high-redshift spirals. But the low metallicities of the damped Lya galaxies cast some doubts on 
this interpretation. A thick disk similar to that of our own Galaxy also appears to be inconsistent 
with the observed abundance distribution in damped Lya galaxies, since the mean metallicity of 
stars in the thick disk is around [Fe/H]=— 0.6 with a range of —1.4 to (cf. Pardi, Ferrini, & 
Matteucci 1995). However, Beers & Sommer-Larsen (1995) suggested that the [Fe/H] distribution 
of thick disk stars probably extends to as low as —2. Since existing studies of thick-disk stars 
arc still limited to within a few kpc of the disk, the known metallicity distribution may be 
biased toward higher values. Thus a thick-disk interpretation of the damped Lya abundances 
remains viable, especially in light of the (somewhat preliminary) evidence that damped Lya 
galaxies exhibit kinematics expected of rotating structures (Wolfe 1995). Detailed accounts of the 
properties of the thick disk of the Milky Way may be found in Majewski (1993). 

On the other hand, the damped Lya galaxies have metallicities resembling those found for 
Galactic halo stars and globular clusters. Coincidentally, the range of [Fe/H] spanned by the 
damped Lya galaxies at 2; ~ 2 — 3 is nearly identical to that spanned by Galactic globular clusters. 
This suggests the interesting possibility that damped Lya absorbers may represent a phase when 
the spheroidal component of high-redshift galaxies is being formed. Such a possibility has, in fact, 
been suggested previously based on different arguments (Lanzctta et al. 1995; Wolfe 1995). The 
Fe-abundance distribution of damped systems is also consistent with that of old stellar populations 
in local group dwarf spheroidal galaxies (cf. Da Costa 1992). Dwarf spheroidals contain old 
to intermediate age stellar populations but with little gas, and the metallicity distributions of 
their stellar populations indicate episodic star formation in their history. Presumably, most or 
all of them contained interstellar gas in their remote past that would be sufficient to produce 
damped Lya absorption lines. The more gas-rich dwarf irregular galaxies (e.g., LMC and SMC) 
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also appear to contain old and sometimes intermediate age stellar populations, in addition to 
young stellar populations indicating ongoing or recent star formation. The metallicities of the 
old stellar populations in these galaxies are also similar to those seen in damped Lya galaxies. 
These considerations suggest that damped Lya galaxies may represent gas-rich dwarf galaxies or 
sub-galactic fragments which are still collecting themselves to form more massive galaxies. More 
discussion about the nature of damped Lya galaxies may be found in §5. 

4.3. Abundance Ratios: Nucleosynthesis vs Dust Depletion 

In figure 23 we show the abundance ratios of Si, S, Zn, Cr, Mn, and Ni to Fe plotted against 
[Fe/H]. Again, except for the [Cr/Fe] panel, the Cr abundance has been used in place of the Fe 
abundance when the latter is not available. This is justified because all the Lya galaxies for which 
we have both Fe and Cr measurements indicate [Cr/H]~[Fe/H]: ([Cr/Fe]) = 0.04 it 0.10. The 
dispersion in [Cr/Fe] is consistent with the typical measurement uncertainty of ~0.1 dex. Most of 
the abundance limits for C, N, O, and Al in Table 16 are not strong enough to constrain the origin 
of these elements, so they will only be considered in special cases. 

In principle, elemental abundance ratios allow one to infer what kind of nucleosynthetic 
processes may be responsible for the enrichment of the interstellar medium. They are important 
for constraining models of stellar nucleosynthesis and the shape of the IMF. This is because 
different elements have different nucleosynthetic origins and different enrichment time scales. For 
example, essentially all the oxygen in the Galaxy and much of the observable a-process elements 
(including Si, S, Ca and Ti) are produced by massive stars (M > IOMq) through explosive nuclear 
burning during Type II SN explosions (cf. Wheeler et al. 1989). The corresponding time scale is 
very short: < 2 x 10'' yrs. In contrast, most of the Fe-peak elements (including V, Cr, Mn, Fe, Co, 
Ni, and possibly the lighter elements Sc and Ti, and the heavier elements Cu and Zn) are believed 
to be produced by Type la SN^ (cf. Wheeler et al. 1989), which evolve from low to intermediate 
mass stars (Mq < M < 8Mq). Thus the time scale for producing the bulk of Fe-group elements 
should be at least as long as the lifetime of the progenitor stars of Type la SN, which is > 10* — 10^ 
yrs, and is possibly much longer depending on the uncertain time scale between the formation of 
the white dwarf and the ignition of SN la explosion. For example, Yoshii, Tsujimoto, & Nomoto 
(1996) deduced an "effective" lifetime of ~ 1.5 Gyrs for Type la SN progenitors based on the 
analysis of Galactic stellar abundance data. The same low-to-intermediate mass stars also make 
much of the solar N (and possibly C) through mass losses on comparable time scales. Hence 
the production of C, N, and Fe-group elements is expected to lag behind that of O and a-group 



^ Using the Type II supernova models of Woosley & Weaver (1995) and a simple chemical evolution model, Timmes, 
Woosley, & Weaver (1995) found that only 1/3 of the solar Fe abundance comes from Type la supernovae. However, 
they also found that better agreement with most of the observed abundance data can be achieved if the Type II 
supernova yields of Woosley & Weaver are systematically reduced by a factor of 2. 
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elements following an episode of star formation. 

However, in the case of elemental abundances derived from gas phase absorption, one also 
needs to worry about the additional effect of dust depletion if a significant amount of dust exists in 
these damped Lya galaxies. Refractory elements such as Si, Mn, Fe, Cr, Ni, etc. are preferentially 
incorporated into dust grains, which depletes their abundances in the gas phase. That different 
elements deplete in different amounts significantly modifies the relative abundances of these 
elements in the gas phase. 

In order to distinguish between the effects of miclcosynthetic processes and dust depletion, we 
compare the abundance ratios observed in damped Lya galaxies with those seen in Galactic disk 
and halo stars, and with those determined for diffuse Galactic ISM clouds. The former serves to 
illustrate the effects of nucleosynthesis at different stages of chemical evolution of the Milky Way, 
while the latter illustrates the effects of dust depletion. Granted, there is no guarantee that any of 
the damped Lya galaxies would have the same chemical enrichment history as the Milky Way, or 
that the dust grains in these high-redshift galaxies would have properties similar to Galactic dust, 
but these comparisons should nevertheless be useful in providing some general guidance. 

4.3.1. Can Dust Depletion Explain the Observed Abundance Ratios? 

Figure 24 shows the abundance ratios in diffuse ISM clouds in both the Galactic disk and halo 
from Table 17. Reviews of this subject may be found in Jenkins (1987) and Savage &; Sembach 
(1996). Since it is believed that the intrinsic elemental abundances of the local interstellar medium 
are solar, the spread in [Fe/H] (the horizontal axis in figure 24) represents different levels of dust 
depletion in the clouds rather than variations in the intrinsic metallicity of the clouds (as is the 
case for figure 23) . The departure of the observed elemental abundances in ISM clouds from the 
solar mixture is generally considered to be caused by condensation of refractory elements into solid 
forms (grains). In the diffuse Galactic ISM, Fe, Cr, and Ni are among the most heavily depleted 
elements; Si and Mn are somewhat less depleted; while Zn, S and C, N, O are nearly unaffected 
by dust depletion (Jenkins 1987). These depletion characteristics show up clearly in figure 24. 
Comparing figure 23 and figure 24, one finds good agreement in the [Cr/Fe] and [Ni/Fe] ratios. 
The [Si/Fe], [S/Fe] and [Zn/Fe] ratios in damped Lya galaxies are consistent with the effects of 
dust depletion (figure 24) in the direction of the deviations from the solar ratios, but the amounts 
of deviations are much less in the damped Lya galaxies. This may be explained if the dust-to-gas 
ratio in damped Lya galaxies is lower than the Galactic vahic or the mean gas density in damped 
Lya galaxies is lower; the latter follows from the fact that the degree of elemental depletion is 
observed to correlate with the mean gas density along the line of sight in the Milky Way (see the 
review in Jenkins 1987). However, the Mn/Fe ratio found for damped Lya galaxies is completely 
opposite to what is expected from the dust depletion effect. Another strong piece of evidence 
against a dust depletion interpretation of the observed abundance ratios in damped Lya galaxies 
is their N/0 ratio (see below). 
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The N I AA1199.55, 1200.22, 1200.71 triplet occurs in the Lya forest so it is generally not 

possible to make reliable measurements of their column densities. However, in two damped Lya 
systems (see Table 16), the clear absence of N I absorption and the presence of saturated O 1 
absorption provide firm upper limits on their [N/0] ratio of < —0.76 and < —0.31 dex, respectively. 
Stronger limits can be found for [N/Si]: < —1.31 and < 1.05, respectively. In addition, we found 
—1.21 <[N/Si]< —0.53 for several components in the z^amp = 2.8110 system toward Q 0528—2505 
(see §3.7). Since both Si and O are believed to come from massive stars, and indeed, are observed 
to have solar abundance ratio in Galactic halo stars (cf. Wheeler et al. 1989) and in metal-poor 
dwarf galaxies (Thuan, Izotov, &: Lipovetsky 1995), it is probably safe to assume that [Si/O]=0 in 
damped Lya galaxies. We then find N/0 at least a factor of 3.3 to 20 times lower than the solar 
ratio in these damped Lya galaxies. 

The abundances of N and O in ISM clouds have been measured toward many sightlines 
using the Copernicus satellite by York et al. (1983), whose data were reanalyzed by Hibbert, 
Dufton, Sz Keenan (1985) and by Keenan, Hibbert, Sz Dufton (1985) using improved /-values. 
Similar determinations have been made for a few clouds using the GHRS on board the HST 
(Cardelli, Savage, & Ebbets 1991). Since these measurements were made using the extremely weak 
intersystem transition lines of N I] A1160 and O I] A1335 rather than the strongly saturated N I 
A1200 triplet and O I A1302 transition, they do not suffer significantly from the usual saturation 
problem associated with the strong lines. Their results indicate that the depletion of N and O in 
diffuse ISM clouds is almost always less than 0.3 dex, and in addition, N is usually less depleted 
than O. Hence, it is extremely difficult to explain the observed N/O ratios in damped Lya galaxies 
with dust depletion if the intrinsic N/O ratio is solar. 

Thus we conclude from the above discussion that dust depletion alone cannot explain the 
abundance ratios observed in damped Lya galaxies. 

4-3.2. Can Nucleosynthesis Explain the Observed Abundance Ratios? 

Figure 25 shows the abundance ratios found in Galactic disk and halo stars (see §4.1.2 for 
references). The overabundance of Si to Fe relative to the solar ratio in halo stars ([Fe/H]< —1), 
[Si/Fe]~ 0.3 — 0.4, is now well established, and is thought to reflect the yield of these elements 
from Type II supernovae. A similar overabundance is seen for S (figure 25) and other a-elements 
(cf. Wheeler et al. 1989). All the other elements displayed in figure 25 (Cr, Mn, Fe, Ni, Zn) are 
Fe-peak elements. The fact that the Fe-peak elements are believed to be formed under nuclear 
statistical equilibrium and thus share a common origin may explain why their abundance ratios in 
the halo stars are so similar to the solar values. The only exception is Mn, with ([Mn/Fe]) ^ —0.3 
for halo stars. This "odd-even" effect, namely, the relative underabundance of odd-Z elements 
(where Z is the atomic number) compared to the abundances of even-Z elements of the same 
nucleosynthetic origin, is fairly well established observationally. One theory of why this occurs is 
that the production of odd-Z elements depends on the neutron excess in the nuclear fuel, which 
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in turn depends on the initial metallicity of the stellar composition (Truran &; Arnett 1971). 
Calculations of explosive nucleosynthesis in metal poor stars by Truran & Arnett indicate that the 
yield of both odd-Z nuclei and the neutron-rich isotopes of even-Z nuclei are significantly reduced 
relative to more mctal-rich stars. 

Wc find remarkable agreement in the Si/Fc, S/Fc, Cr/Fc, and Mn/Cr ratios between damped 
Lya galaxies and those seen in Galactic stars in the same metallicity range. The difference in 
Ni/Fe between figures 23 and 25 may not be real and may have to do with the uncertain /-values 
of the Ni II transitions (cf. Morton 1991). Since the stellar abundances are not determined from 
the same Ni II lines (the stellar abundances are from Ni I lines) as for the damped Lya galaxies, 
a systematic error in the Ni II /-values could explain the nearly constant offset of [Ni/Fe] from 
the zero point in damped Lya galaxies. This may also explain why the Ni/Fe ratios in damped 
Lya galaxies more closely resemble those seen in diffuse ISM clouds (figure 24) as in both cases 
the abundances of Ni are derived from the same Ni II transitions with the same /-values. We 
consider the Ni/Fe ratio unsuitable for the purpose of discriminating a nucleosynthetic origin of 
the damped Lya abundances from a dust depletion origin. More reliable determinations of the 
Ni II /-values in the future would make the Ni measurements more useful. We also note that, 
although in most damped Lya systems the Al II A1670 absorption lines are saturated, the Al II 
absorption in the z^amp = 2.8443 system toward Q 1946+7658 and in the Zdamp = 4.0803 system 
toward Q 2237—0608 appears weak enough that the estimated A''(A1 II) 's probably do not suffer 
from significant saturation effects. The ratios, [Al/Si]= —0.30 and —0.35 (Tables 11 and 15), are 
consistent with observations of Galactic halo stars (cf. Wheeler et al. 1989). Incidentally, Al 
and Si are both produced during hydrostatic carbon and neon burning in massive stars, but Al 
is an odd Z element. Accordingly, the underabundance of Al relative to Si is consistent with the 
odd-even effect. 

The only significant difference between figure 23 and figure 25 is in the Zn/Fe ratio. A Zn/Cr 
overabundance, or equivalently, a Zn/Fe overabundance (since [Cr/Fe]~0 in damped Lya galaxies) 
was known since the first studies of damped Lya abundances (Meyer, Welty, & York 1989; Pettini, 
Boksenberg, & Hunstead 1990), and was interpretated to indicate that some Cr (and Fe) is 
depleted by dust grains from the gas phase. We will have more discussion of this in §4.3.3. 

In general, we consider the good agreement in the Si/Fe, S/Fe, Cr/Fe, and Mn/Fe ratios 
between damped Lya galaxies and Galactic halo stars strong evidence for a Type II SN enrichment 
origin of these elements without significant modifications by dust grains. In particular, the 
observed Mn/Fe ratios require a Type II SN-cnrichment explanation since they cannot be explained 
otherwise. Another strong argument for SN II enrichment in damped Lya galaxies is the N/0 
ratios (see below). 

In §4.3.1, we argued that the observed low N/0 ratios in several damped Lya galaxies are 
inconsistent with the interpretation of dust depletion. On the other hand, these low N/0 ratios 
are very similar to those found in Galactic halo stars. There was some debate about whether 
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[N/Fe] in halo stars is more close to zero or less than zero (see discussion in Wheeler et al. 1989). 
But since [0/Fe] is known to be about +0.4 dex in halo stars, it is safe to say that [N/0]< —0.4 
dex at [Fe/H]< —1. Similarly, observations of metal-poor ([0/H]< —1) external galaxies indicate 
[N/0]< -0.6 (Pagel et al. 1992; Skillman k, Kennicutt 1993; Thuan et al. 1995). Such low 
N/0 ratios are commonly interpreted in terms of the different nucleosynthetic origin of these 
elements. While O is almost exclusively produced by short-lived massive stars, N is thought to 
have two sources. For pristine gas or gas with very low metallicity, primary N production is 
made by intermediate mass stars (3Mq < M < 8©) during the third drcdgc-up episodes, which 
bring carbon-rich material from the He-burning shell into the hydrogen-burning shell (Rcnzini & 
Voli 1981). The N produced this way is primary because it does not depend on the initial metal 
content of the gas. The long evolution time scale of the intermediate mass stars necessarily means 
that there should be a time delay of several hundreds of million years between the injection of O 
into the ISM and that of primary N. Secondary production of N can occur in stars of any mass 
during main-sequence burning, when the initial metallicity in the star is high enough to provide 
the seed C and O nuclei. Although the above popular interpretation is challenged by some recent 
observational results (cf. Thuan et al. 1995), which suggest that significant amounts of primary N 
may be produced by massive stars, it remains true that the observed N/0 ratio in halo stars and in 
metal-poor external galaxies is low and comparable to that observed in damped Lya galaxies. We 
thus believe that the low N/0 ratios (and the low Mn/Fe ratios) found for damped Lya galaxies 
are indications that massive stars are the primary pollutants of the ISM in these galaxies. This 
probably happens because star formation in these galaxies has occurred shortly before they were 
observed, so that only massive stars have had time to evolve to Type II SN. Given this time scale, 
we expect to see a-elements overabundance relative to Fe-peak elements. This is clearly seen in 
figure 23. In fact, the observed overabundances of Si to Fe and S to Fe, and the under abundance 
of Al to Si and Mn to Fe can be entirely explained away with the above nucleosynthesis argument. 
The only remaining question is: is there a need for dust depletion on top of the nucleosynthesis 
pattern? 



4-3.3. Dust Depletions On Top of Nucleosynthesis Pattern? 

The preceding discussion demonstrates that dust depletion alone cannot explain the observed 
elemental abundances in damped Lya galaxies. Rather, with the exception of Zn, all the observed 
abundances are nicely explained if Type II supernovae are the only significant sources of pollution 
in these galaxies. Could the Zn/Fe ratios, then, be suggesting that there is some dust depletion 
effect on top of the nucleosynthesis pattern? 

In Galactic stars the Zn/Fe ratio is found to be solar at all metallicities (see figure 25). On 
the other hand, the observed Zn/Fe ratios in damped Lya galaxies are essentially all above the 
solar ratio. Since Zn is almost unaffected by dust in warm diffuse ISM clouds, while Fe is among 
the most heavily depleted elements, a modest amount of Fe depletion would explain the observed 
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Zn/Fe ratios in damped Lya systems. This was indeed the argument that was put forward 
previously to explain the observed Zn/Cr (equivalent to Zn/Fe) ratios (cf. Pettini et al. 1994). 
However, this interpretation does not appear to be consistent with the other observed ratios. For 
the damped Lya systems in which both Zn and Fe abundances are measured, ([Zn/Fe]) = +0.39. 
In other words, if the above dust-depletion explanation for the observed Zn/Fe ratios is correct, 
the gas phase Fe abundance in damped Lya galaxies must have been on average depleted by 
about 0.4 dex (or more if Zn is also depleted). To explain the other ratios shown in figure 23, 
this requires Si, S, Cr, Mn being depleted by the same amount in order to maintain the expected 
Si/Fe, S/Fe, Cr/Fe, and Mn/Fe ratios from halo star observations and from theories of Type II SN 
enrichment (recall that the observed Mn/Fe and N/0 ratios in damped Lya galaxies require that 
the enrichment process be dominated by SN II). This is quite difficult to understand given the 
observed depletion patterns in Galactic ISM clouds. As discussed in Sembach & Savage (1996), 
even in warm diffuse halo clouds for which elemental depletion is the least, there is a differential 
depletion factor of 0.44 dex between Si and Fe (i.e., Fe is more depleted than Si by 0.44 dex). The 
difference in the Si and Fe depletion widens for cooler, denser clouds. The problem is especially 
acute for S since in diffuse ISM clouds S is essentially undepleted (see Table 17). Given that Fe is 
among the most heavily depleted elements in the Galactic ISM, it is difficult to understand how S 
and Fe could have the same depletion, unless the depletion property of dust in the damped Lya 
galaxies is drastically different from Galactic dust. However, if the latter is true, then the basis for 
interpretating the Zn/Fe ratio to indicate dust depletion may no longer be valid. 

The degree to which an element is depleted by dust grains in the Galactic ISM is known to 
roughly anticorrelate with its condensation temperature (Field 1974; Jenkins 1987), which is the 
temperature at which half of the atoms condenses into solid forms of some sort under thermal 
equilibrium conditions. To further look into the issue of dust depletion, we plot the abundances 
of the elements found in damped Lya galaxies as a function of their condensation temperature 
using the data given in Table 16. Rather than plotting the abundance of individual elements for 
each damped Lya galaxy, which would make the figure incomprehensible because of the large 
scatter, we chose to plot the elemental abundances relative to the Fe abundance in the damped 
Lya galaxies in order to sec the mean trend. The result, shown in figure 26, is both quantitatively 
and qualitatively different from the same relation seen in ISM clouds (cf. Jenkins 1987; Savage & 
Sembach 1996). The most notable difference is in the behaviour of Mn and N. The implication is 
that there is no significant evidence for dust depletion to be the origin of the observed abundance 
patterns in damped Lya galaxies. This is in agreement with similar findings for two individual 
damped Lya galaxies previously studied by Wolfe (1995) and by Prochaska & Wolfe (1996). 

Another way to look for the effect of dust depletion is to examine the variations of elemental 
abundance ratios among the different velocity components in a given system. It is well established 
that the degree of dust depletion in ISM clouds generally varies with their locations (environment) 
in the Galaxy, with cool, diffuse disk clouds showing much larger depletion than the warm, 
diffuse halo clouds (cf. Sembach Sz Savage 1996). Such a dependence of depletion on environment 
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probably has to do with the modification or destruction of dust grains by supernova shocks. If the 

damped Lya galaxies contain dust grains that are affected by similar energetic processes, then 
one might expect to see large variations in the elemental abundance ratios from one absorbing 
cloud to another. The best elements for this test are probably (S, Zn) and (Cr, Fe, Ni) since the 
two groups of elements are at the two extremes of dust depletion and because they are generally 
measurable in damped Lya systems. In practice, however, such a test is difficult to carry out since 
one must rely on the weak lines to derive accurate column densities, which require high resolution, 
high S/N measurements. Prochaska & Wolfe (1996) examined the relative abundances of Si, Fe, 
Ni, Cr, and Zn in several components of the Zdamp = 2.462 damped Lya absorption system toward 
Q 0201+3634 and concluded that there is little dust present in the absorbing gas. Among all the 
damped Lya systems studied in §3, we consider the Zdamp = 2.8110 system toward Q 0528—2505 
the best suited for this purpose since the S II A1250, 1253 lines are well measured and we have 
good measurements for six Ni II lines (see Table 8). Equally important is the fact that these ions 
are well measured across a ~ 300 km s~^ interval encompassing many components. The three 
panels to the left in figure 27 show the profiles of Na{v), the column density per unit velocity 
interval, for S II and Ni II derived from the average of two S II lines and of six Ni II lines in the 
system. The profile for Si II A1808 is also shown for comparison. The three panels to the right in 
figure 27 shows the Na{v) ratios of these ions. Amazingly, the S Il/Si II, S II/Ni II, and Si II/Ni 
II ratios are nearly constant over the Av ~ 300 km/s interval where these ratios are relatively well 
measured. Given that these ratios can be affected by many factors, including intrinsic relative 
abundance variations, ionization effects, as well as dust depletion, we find the relative constancy of 
these ratios quite remarkable. In comparison, the S II/Ni II ratio varies from ~ 6 in warm diffuse 
halo clouds to > 200 in cool disk clouds in the Milky Way (see Sembach & Savage 1996). The 
above result becomes even more remarkable in light of the fact that the Zdamp = 2.8110 system 
toward Q 0528—2505 is the only damped Lya system known to contain molecular hydrogen (Foltz 
et al. 1988; Songaila &; Cowic 1996), which makes it more likely to contain dust. However, in this 
damped Lya system the fraction of molecular hydrogen relative to neutral hydrogen is well below 
the Galactic value for sightlines with similar values of A''(H I). 

In conclusion, wc find it difficult to interpret the observed Zn/Fe ratios (or Zn/Cr ratios) in 
damped Lya galaxies as being a consequence of dust depiction. 

4-3.4- Concluding Remarks 

To summarize the preceding discussion, we favor the interpretation that the abundance ratios 
in damped Lya galaxies are mostly or entirely a consequence of Type II SN nucleosynthesis 
without significant modification by dust depletion. The evidence for this conclusion is the 

much-lower-than-solar N/0 ratio, the a-clcment overabundance relative to Fe peak elements, and 
the undcrabundance of Mn relative to Fe and Al relative to Si (odd-even effect). Such abundance 
patterns are, in some sense, expected in the early stages of galactic chemical evolution based on 
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what we know about the Milky Way (cf. Timmes, Lauroesch, &; Truran 1995), and were indeed 
noted previously in individual damped Lya systems (e.g., Lu et al. 1995b: Meyer, Lanzetta, & 
Wolfe 1995; Steidel ct al. 1995; Wolfe 1995; Pcttini et al. 1995a; Lu ct al. 1996a; Prochaska & 
Wolfe 1996). The present study, however, puts the above conclusion on much firmer ground. 

The Zn/Fc overabundance found for damped Lya galaxies can be explained without invoking 
dust if there is an additional source of Zn (particularly ^^Zn, see figure 5 of Timmes et al. 1995) 
in SN II besides that synthesized under nuclear statistical equilibrium. Such an alternative 
production mechanism for Zn might occur in the neutrino-driven wind discussed by Hoffman et 
al. (1996). One difficulty with this interpretation is that a Zn/Fe overabundance is not observed 
during the early history of the Milky Way (figure 25). In that regard, it would be extremely 
valuable to obtain the Zn/Fe ratio for some metal-poor stars in external galaxies in order to see 
if Zn/Fe is universally solar. It is perhaps significant that we do not yet have a clear theoretical 
understanding of why Zn should track Fe so well in Galactic disk and halo stars. For example, 
using the Type II SN yield calculations of Woosley Sz Weaver (1995), both Timmes et al. (1995) 
and Malaney & Chaboyer (1996) found a much lower-than-solar Zn/Fe ratio (by a factor of 
3-7) for plausible assumptions about the stellar initial mass function (sec Arnett 1995 for a 
more in-depth discussion). If indeed the observed Zn/Fe ratios in damped systems are intrinsic 
to the nucleosynthetic processes, it will have significant implications for the theory of stellar 
nucleosynthesis. At this point, we consider Zn an anomalous element. It is important to note 
that, even if the observed Zn/Fe ratios are due to dust depletion (which we consider unlikely), 
the amount of Fc depletion should be no higher than 0.4 dex on average, assuming that Zn is 
undepleted. In other words, the amount of Fe locked up in grains should be less than 60% of the 
total. This is in contrast to the depletion of Fe in ISM clouds in which up to more than 99% of 
the Fe may be locked up in grains. It should also be noted that a 0.4 dex depletion of Fe would 
not explain all the differences between the metallicity distribution of damped Lya galaxies and 
that of the Milky Way disk stars (figure 22). 

Regardless of whether or not dust is present in the damped Lya galaxies, the observed N/0 
ratios and Mn/Fe ratios strongly indicate a nucleosynthetic origin. This suggests that the chemical 

enrichment process in these damped Lya galaxies must have occurred so recently prior to when 
they were observed that only Type II SN have made significant contributions to the pollution of 
their interstellar material; low to intermediate mass stars could not have had enough time to evolve 
and to make their presence known by dumping their nucleosynthetic products into the interstellar 
space. Hence the star formation process should not have proceeded more than several hundred 
million years from the time corresponding to their redshifts. One way to avoid this conclusion is if 
the stellar initial mass function in these young galaxies had a low mass cutoff so that there simply 
aren't that many low mass stars to make Type la SN, even if they had enough time to do so. 
Observations of elemental abundance ratios in lower redshift and more metal-rich damped Lya 
galaxies should help to clarify this issue, since we expect SN la to catch up eventually. 



-30- 



4.4. Discussion 

The issue of dust is critical for a complete and accurate understanding of damped Lya 
galaxies. Fall and collaborators (Fall & Pei 1989; Fall, Pei, & McMahon 1989; Pei, Fall, & Bechtold 
1991) found marginally significant evidence that quasars with known damped Lya galaxies in 
the foreground are statistically redder than those without. This was taken as evidence for the 
presence of a modest amount of dust in damped Lya galaxies. The dust-to-gas ratio inferred 
for the damped Lya galaxies is roughly 1/10 of that in the Galaxy. However, the observational 
evidence for this reddening has been challenged recently (Foltz 1996, private communication; 
Storrie-Lombardi 1996, private communication). The only other independent evidence for the 
presence of dust in damped Lya galaxies was thought to be the Zn/Cr or Zn/Fe ratios (cf. Pettini 
et al. 1994). But as we argued earlier, this interpretation has problems explaining many of the 
other observed abundance ratios. 

Fall &; Pei (1995; also Pei &; Fall 1995) argued that most of the more dust-rich galaxies 
have been missed in existing surveys. This is because previous surveys for damped Lya galaxies 
were based on quasars selected optically, hence quasars with dusty foreground galaxies would be 
preferentially missed owing to dust obscuration. Fall & Pei (1995) estimated that up to 70% of 
the bright quasars at z ~ 3 could be missing from optically selected samples. In principle, such a 
selection effect could occur even if there is no significant amount of dust in the existing sample of 
damped Lya galaxies. If this is the case, then the damped Lya galaxies studied here will not be 
representative of all high-redshift galaxies. In particular, the mean metallicity as given by [Fe/H] 
would be biased towards lower values, making them appearing (on average) more metal poor than 
they really are. The abundance ratios may also be affected. A critical test of this hypothesis is to 
conduct a new survey for damped Lya galaxies using a sample of high redshift quasars selected 
from a complete sample of radio sources regardless of their optical brightness. Such a program is 
being pursued. 

5. OTHER PROPERTIES AND THE NATURE OF DAMPED LYMAN- ALPHA 

GALAXIES 

5.1. Electron Densities 

The column density ratio of C II*/C II is a density indicator (Bahcall Sz Wolfe 1968). The 
lower level of the C II* A 1335 transition is a fine structure level of the ground state of the 
C II atom. Under normal interstellar conditions, the population balance between the ground 
level and the fine structure level of the C II ground state is determined by collisional excitation 
with electrons and by spontaneous radiative decay. Table 18 summarizes the C II* and C II 
measurements in damped Lya systems studied here and in Lu et al. (1996a). Upper limits on the 
electron density are derived using the formula and coefficients adopted in Tripp et al. (1996) by 
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assuming a temperature of 5000 K. The derived electron densities are upper limits because only 
lower limits to A^(C II) are available in all systems considered, and because other mechanisms 
might also contribute significantly to the excitation (e.g., excitation by microwave background 
photons; see §6). These limits are consistent with the mean electron density of ~ 0.07 cm~^ 
estimated for the Galactic ISM (cf. Reynolds 1991; Spitzer & Fitzpatrick 1993). For reference, 
the derived upper limits are on average factors of 2-3 lower if we estimate N{C II) by assuming a 
solar C/Si ratio. 

5.2. Ionization 

The relative strength and velocity distribution of absorption from ion species with different 

ionization properties can provide insight into the ionization structure and the spatial distribution 
of the absorbing gas. Tabic 19 summarizes such information for some ions of interest based on 
data presented in this paper and other similar studies with Keck HIRES (see references in the 
table). We make the following comments: 

(1) In all cases where such information is available, the Al III absorption profiles always 
have similar appearance as the low ionization lines. Additionally, the Al II lines are always much 
stronger than the Al III lines, indicating N{A\ II)>> N{Al III). These results suggest that Al III 
is likely to be produced in the same physical region as the low ion gas, and that most of the gas is 
neutral rather than ionized (as expected from the large H I column densities of the systems). This 
could happen if Al III comes from an ionized shell surrounding the neutral gas which produces the 
bulk of low ions. 

(2) The Si IV and C IV absorption line profiles almost always resemble each other, and they 
usually have a very diflcerent appearance from the low ion absorption lines. This suggests that the 
bulk of the high ions are probably produced in physical regions distinct from that containing the 
low ions. The high ions could arise from low-density ionized gas (halo clouds?) surrounding the 
main structure which gives rise to the low ionization species. 

The profile differences between low-ion and high-ion absorption lines have been noted 
previously in individual damped Lya systems (eg. Wolfe et al. 1994; Lu et al. 1995b). The 

current discussion makes these trends much clearer. The results are reassuring in that, in general, 
the metal abundances derived from low ionization species without corrections for ionization eff'ects 
should be substantially correct. This is consistent with theoretical expectations (Viegas 1995) 

5.3. Kinematics 

The profiles of metal absorption line in damped Lya systems contain information about the 
kinematics of the gas. The total velocity spread of the absorption lines may allow one to estimate 
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the mass of the objects if they are caused by random motions of the absorbing clouds confined by 

gravitational potentials. In practice, the velocity spread may be affected by ordered motions (e.g., 
rotation, infall/outflow) and by energetic events such as supernova explosions. One needs also 
to take into account the fact that often the quasar sightlines do not probe the full gravitational 
potential of the absorbing objects, depending on the impact parameters. These considerations 
make it difficult to infer the mass of the absorbing galaxies. 

On the other hand, analyses of the profiles of the metal absorption lines have yielded a 
potentially very exciting result. Wolfe (1995; see also Lanzetta & Bowen 1992) found preliminary 
evidence that the weak low-ionization metal absorption lines in damped Lya systems often show 
an edge-leading asymmetry that is characteristic of absorption from rotating gaseous disks. It 
should be pointed out that only weak, unsaturated low-ionization absorption lines can be used 
to trace reliably the velocity structure of the neutral gas in galactic disks. Strong, saturated 
absorption lines are not appropriate for this purpose for two reasons: (1) one loses the information 
on the velocity structure of the absorption in the saturated part; (2) strong absorption lines are 
more sensitive to trace amounts of gas so they are more easily contaminated by diffuse clouds 
that are not associated with galactic disks (e.g., halo clouds). The implication of Wolfe's finding 
is that damped Lya galaxies may be high-rcdshift spiral disks. If this conclusion is confirmed 
with a much larger data set, it will not only provide strong evidence for the physical nature of the 
damped Lya galaxies, but will also become one of the most significant constraints on cosmological 
models of structure formation in the early universe. Data contained in this paper are currently 
being analyzed by Wolfe and collaborators for the purpose discussed above. 



5.4. Star Formation and Stellar Initial Mass Function 

The rate and form of star formation and the initial mass function (IMF) of the stars formed 
are fundamental parameters for understanding galactic chemical evolution. The failure to detect 
Lya emission from damped Lya galaxies may suggest that damped Lya galaxies generally do not 
have high star formation rates. For example, Lowenthal et al. (1995) surveyed seven quasar fields 
with known foreground damped Lya absorbers and failed to detect any of the galaxies in Lya 
emission down to quite stringent limits. The implied star formation rates are < IMq yr~^ for 
normal IMF, assuming negligible extinction by dust. Of course, if a significant amount of dust is 
present in these galaxies, the true star formation rate could be much higher. However, the Hu et 
al. (1993) observations of Ha emission from three damped Lya galaxies limit the star formation 
rates to less than 20Mq yr~^ in those galaxies. 

On the other hand, the lack of Lya emission from a damped Lya galaxy does not necessarily 
mean low star-formation activity, even if there is no significant extinction by dust. As demonstrated 
by Valls-Gabaud (1993), the equivalent width of Lya emission from a galaxy is a strong function 
of the age of the stellar population. In the case of a star burst, the strength of the Lya emission 
decreases rapidly after the burst, essentially because the massive stars (O and early B stars) 
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responsible for the Lya emission are short-lived. His model calculations suggest that the Lya 
photons disappear from the galaxy in ~ 10^ years after the burst and the galaxy spectrum will 
even show stellar absorption in Lya at later times. The case of continuous star formation is 
qualitatively different, with the galaxy showing relatively strong Lya emission (equivalent width 
> 100 A) at all times because O, B stars are continuously being formed. If the lack of strong 
Lya emission in damped Lya galaxies is indicative of the bursting nature of their star formation 
process, it may have significant implications for their nature. There is now strong evidence that 
normal spiral disks tend to have relatively constant star formation during their lifetimes while 
dwarf galaxies (dwarf irregulars and dwarf sphcroidals) tend to have bursts of star formation (cf. 
Kennicutt 1995). A bursting nature of star formation for damped Lya galaxies, if established, 
would favor them being dwarf galaxies or subgalactic structures. 

The form of star formation in damped Lya galaxies may also be inferred from the elemental 
abundance ratios. In particular, the relative abundance of O and a-process elements to those 
of the Fe-group elements is fairly sensitive to the intervals between episodes of star formation, 
essentially because of the ^ 1 Gyr delay between the almost-instantaneous enrichment from SN 

11 and that from SN la. If some of the damped Lya galaxies in our sample began with a strong 
burst of star formation which was followed by a long quiescent period, one would expect its Si/Fe 
ratio to drop below the SN II expectation (i.e., ~2-3 times the solar value) after Ai >> 1 Gyrs 
when Type la SN start to inject Fe into the ISM. Thus a uniform distribution of Si/Fe ratios 
may well indicate that the star formation process in these galaxies is relatively uniform. The 
existing data (figure 23) docs show a fairly uniform Si/Fe ratio: ([Si/Fe]) = 0.36 it 0.11 based on 

12 measurements (note that the scatter of 0.11 dex is comparable to the typical measurement 
uncertainty of 0.1 dex). However, it is premature to make any firm conclusions because of the 
relatively small number of data points, and because we don't know when each galaxy was formed. 
The age span of the existing sample galaxies with Si/Fe measurements, ~ 3 Gyrs, is probably 
too small to even make a statistical inference about when these galaxies were formed. If the 
uniformity of the Si/Fc ratio persists to much lower redshifts in a much larger sample, then it may 
seriously constrain the history of star formation in these galaxies since at least some of the very 
low-redshift galaxies should have formed much longer than 1 Gyrs ago. A smooth star formation 
process would be more typical of disk galaxies than dwarfs. Alternatively, a constant Si/Fe ratio 
may be maintained by a top-heavy IMF. Dust depletion, which is expected to become significant 
at low enough redshift, may further complicate the issue. 

The relative constancy in the Si/Fe ratio among the damped Lya galaxies in the sample also 
suggests that the IMF integrated stellar yields of these elements are fairly uniform. If the relative 
yield of a-elements and Fc from Type II SN varies strongly with the progenitor mass, as was 
suggested in some theoretical calculations (e.g., Thielemann, Nomoto, Sz Hashimoto 1996), the 
scatter in the measured Si/Fe ratio may be used to put significant constraints on the variations 
of the stellar IMF in these galaxies over the range of 10-40 Mq (responsible for SN II). We did 
not attempt a quantitative calculation of the IMF constraint because all existing models of SN 
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II suffer from the uncertain nature of the explosion mechanism, which significantly affects the Fe 
yield. 

5.5. Are Damped Lya Galaxies Spirals? 

There is direct evidence that at least some damped Lya absorption systems may be associated 
with large, disk-like structures. Briggs et al. (1989) found that the 21-cm absorption from 
the 

^damp — 2.04 damped Lya galaxy toward PKS 0458—020 shows striking similarity in the 
measurements made with the Arecibo single dish and with the VLBI. This requires the absorbing 
galaxy to cover the entire radio source, which is an extended structure of about 10 kpc size in 
the absorber's reference frame. Kinematic arguments further suggest that the actual size of the 
absorber may be several times bigger. More recently, Djorgovski et al. (1996) identified a galaxy 
at z = 3.150 responsible for a damped Lya absorption system in the spectrum of a background 
quasar. The projected separation between the detected galaxy and the quasar then indicates that 
the absorbing galaxy must be at least as large as 17 kpc, comparable to the sizes of normal galaxy 
disks. 

Another evidence for high-redshift damped Lya galaxies to be spiral disks comes from the 
work of Wolfe (1995), who presented evidence that the absorption profiles of the metal absorption 
lines in damped Lya galaxies generally show an edge-leading asymmetry that is characteristic of 

absorption from rotating gaseous disks. If this result is confirmed with more detailed modeling 
using a large sample, it will constitute the strongest evidence that damped Lya galaxies are spiral 
disks or proto-disks. 

Evidence regarding the nature of damped Lya galaxies also comes from direct imaging 
identification of their low-redshift counterparts. Steidel et al. (1994) presented ground-based 
images of two quasar fields, 3C 286 (Q 1328-1-3045) and PKS 1229—021, containing foreground 
low-redshift {z < 1) damped Lya absorbers, and identified candidate damped Lya absorbing 
galaxies (i.e., no spectroscopic redshifts were obtained). The galaxy responsible for the damped 
Lya absorption toward 3C 286 appears to be a low surface brightness galaxy, while the one 
responsible for the damped Lya absorption toward PKS 1229—021 appears to be a normal galaxy 
somewhat fainter than L*. In another study (Cohen et al. 1996), the damped Lya absorber at 
Zdamp = 0.437 toward the quasar 3C 196 was identified with a spiral galaxy of luminosity roughly 
L* in an HST image, although no spectroscopic redshift confirmation was available. 

More broadly, Steidel and collaborators (see Steidel 1995 and references therein; see also 
Bergeron and Boisse 1991) have been very successful in identifying galaxies responsible for the 
Mg II absorption systems seen in quasar spectra at Zabs < 1-6. The colors and luminosities of the 
identified galaxies are consistent with them being normal Hubble-sequence galaxies, spanning the 
range of E to Irr but mostly early to late type spirals. The implication is that normal galaxies are 
already in place by a redshift of 1.6 with properties not that different from the local population. 
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This is consistent with results from deep surveys of field galaxies (cf. Driver et al. 1995). 

Direct searches for high-redshift galaxies are also starting to produce interesting results, some 
relevant to the nature of damped Lya galaxies. Identifications of z > 3 galaxies in quasar fields 
based on the presence of a redshifted Lyman continuum break (Steidel et al. 1996a,b; Giavalisco 
et al. 1996) have produced a population of what appears to the progenitors of normal, massive 
galaxies that are forming stars of their spheroidal component, although this interpretation is not 
unique. One such "Lyman break" galaxy is confirmed to be a damped Lya galaxy (Steidel et al. 
1996a; Djorgovski et al. 1996), and another is a candidate for a damped Lya galaxy (Giavalisco 
et al. 1996). 

These results suggest that it may not be unreasonable to think that damped Lya galaxies are 
high-redshift spirals or their precursors. In fact, if one takes the age estimates of the disk stars in 
the sample of Edvardsson et al. (1993) literally^ (see figure 23), then it is clear that the Milky 
Way disk was already in place 10-15 Gyrs ago or at z = 2 — 4 for the cosmology adopted here 
(go = 0.1, Hq = 50). Unless our Milky Way is among the first galaxies formed in the universe, one 
may be forced to conclude that (some) spiral disks must exist by the redshift of 2-4 in some form. 
If we insert a sightline through such a galaxy, a damped Lya absorption will clearly result. The 
key question is: what dominates the absorption cross- section at N{m)=W^^ - 10^^ cm at the 
relevant redshifts, spirals or subgalactic structures? 

If the damped Lya galaxies in our sample are indeed disks or proto-disks, we must explain 
why their Fe metallicities are so much lower than the Milky Way disk in its past. One possibility 
is that this is a selection bias due to dust obscuration (Pei &: Fall 1995; see discussion in §4.4). In 
principle, disks or proto-disk could be present at z > 2 without being included in our sample if the 
existing sample of damped Lya galaxies is dominated by a large population of smaller structures. 
Such subgalactic structures could be the building blocks of large, normal galaxies that are still 
falling into each other, or they might survive to the present epoch. If the damped Lya galaxies 
at z = 2 — 4 are predominantly sub-galactic structures, then the number density of damped Lya 
galaxies, dN/dz = 0.3 at z ~ 3 (Lanzetta et al. 1991; Wolfe et al. 1996), requires a space density 
of 322/i5o/i?o Mpc~^, where Rq is the characteristic size (radius in units of kpc) of the subgalactic 
structures with A^(II I)> 2 x 10'^'' cm~^ for a spherical geometry. This may provide significant 
constraints on cosmological models of structure formation. 



Private communications with Dr. Edvardsson suggest that, due to observational uncertainties and uncertainties 
in stellar evolution models, the oldest disk stars in the Edvardsson et al. (1993) sample are not necessarily as old 
as were suggested in the study. Some of the oldest stars in the sample may also be halo stars that happen to have 
orbits very close to the Galactic plane. A detailed discussion of the age of the Galactic disk is beyond the scope of 
this paper. Interested readers are referred to the nice discussion by Majewski (1993). 
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6. TEMPERATURE OF THE COSMIC MICROWAVE BACKGROUND 

Big Bang cosmology predicts that the temperature Tcmb of the cosmic microwave background 
(CMB) radiation should behave like Tcmb oc (1 + As suggested by Bahcall & Wolf (1968), this 
relation is testable since one can estimate the CMB temperature from the relative populations 
of ground-state finc-structTirc levels of suitably selected atoms in high-redshift quasar absorption 
systems. This technique has previously been applied to the C I atom (Meyer et al. 1986; Songaila 
et al. 1994a) and to the C II ion (Songaila ct al. 1994b; Lu et al. 1996a). In §3, we have estimated 
column densities (or limits) of C II and C II* for several damped Lyct absorption systems. Since 
in most cases we only have lower limits to A^(C II) and/or upper limits to N{C II*), the resulting 
upper limits on N{C II*)/N{C II) yield strict upper limits on the temperature of the CMB since 
other mechanisms may contribute significantly to the excitation of the C II ion (e.g., coUisional 
excitations by electrons or hydrogen atoms, direct UV pumping, etc.). 

Table 20 summarizes the available measurements of the temperature of the cosmic microwave 
background at high redshifts, including four based on measurements made in §3, one from Lu et 
al. (1996a), and three from Songaila et al. (1994a, b). In the case of Q 1331+1704 (Songaila et al. 
1994a), both the C I and C I* absorption are detected so it was possible to actually measure the 
excitation temperature rather than providing upper limits. These results are displayed in figure 
28, where the dotted line is the prediction of Tcmb = 2.73(1 + z) from Big Bang cosmology. All 
the existing measurements or upper limits are consistent with the Big Bang prediction. 

We note that, other than the two measurements aX z = 1.776 (which are based on C I*/C I), 
the two upper limits at 2; > 4 are particularly tight. This may be partially attributed to the good 
S/N of the data, but is also due to the fact that at 2; > 4, the peak of the microwave background 
radiation is shifted to a wavelength region very close to the energy separation (158 fim) between 
the fine structure levels of the C II ion ground state. Thus the C II ion is a particularly good 
probe of the microwave background temperature at very high redshifts. It is also significant that 
one of the measurements at z = 1.7660 and the upper limit at z = 4.0803 are very close to the Big 
Bang predictions, suggesting that contributions from other excitation mechanisms are negligible 
in these systems. 

7. SUMMARY 

We have obtained high quality spectra of quasars using the High Resolution Spectrometer 
(HIRES) on the 10m Keck telescope in order to study the elemental abundances in 14 high-rcdshift 
(0.9 < z < 4.0) damped Lya galaxies. The HIRES spectra have a resolution of FWHM=6.6 
km s"^, with a typical S/N~ 40 per resolution element. The chemical elements being studied 
include C, N, O, Al, Si, S, Cr, Mn, Fe, Ni, and Zn, although not every element listed is measured 
for every damped Lya galaxy. In order to ensure accuracy, we measure column densities and 
abundances using only those absorption lines with > 4cr detection significance that are believed 
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to be unsaturated. This allows us to circumvent the difficulty associated with measuring reliable 
column densities using saturated absorption lines. In all other cases, we derive abundance limits 
using either the 4a upper limits (for weaker or undetected lines) or the lower limits (for saturated 
lines) on their column densities. We also use the improved oscillator strengths (/-values) for a 
number of important transitions as given in the compilation of Tripp et al. (1996). These new 
abundance measurements are combined with similar measurements in the literature that are 
believed to be largely free of biases resulting from line saturations in order to investigate the 
chemical evolution of damped Lya galaxies. Where appropriate, we have corrected the earlier 
measurements with the set of improved /-values (Tripp et al. 1996) so that all the measurements 
will be on the same footing. The total sample contains 23 damped Lya absorption systems in the 
redshift range 0.7 < z < 4.4. Our main results are the following. 

1. The damped Lya galaxies have [Fe/H] in the range of —2.5 to —1, or (Fe/H) between 1/10 
to 1/300 solar, clearly indicating that these are young galaxies in the early stages of chemical 
evolution. The 7V(H I)-weighted mean metallicity of the damped Lya galaxies with 2 < z < 3 is 
[Fe/H]=— 1.56, which may be regarded as the cosmic metallicity at {z) = 2.5. There is a large 
scatter (about a factor of 30) in (Fe/H) at z < 3, which we argue is real and which probably 
results from the different formation histories of the absorbing galaxies or a mix of galaxy types. 

2. Comparisons of the distribution of [Fe/H] vs redshift for the sample of damped Lya 
galaxies with the similar relation for the Milky Way disk (i.e., the Galactic age- metallicity relation) 
as defined by the stars in the Edvardsson et al. (1993) sample indicate that the damped Lya 
galaxies are much less metal-enriched than the Galactic disk in its past. Since there is evidence 
from our analyses that depletion of Fe by dust grains in the damped Lya galaxies is unimportant, 
the difference in the enrichment level between the sample of damped Lya galaxies and the Milky 
Way disk suggests that the damped Lya galaxies are probably not high-redshift spiral disks in the 
traditional sense. Rather, they could represent a thick disk phase of galaxies, or more likely the 
spheroidal component of galaxies, or dwarf galaxies. 

3. There is evidence that the mean metallicity in damped Lya galaxies rises sharply at 

z < 3, possibly indicating the onset of significant star formation and chemical enrichment in early 
galaxies. Such an epoch may be identified as the epoch of galaxy formation in the sense that 
galaxies are beginning to form the bulk of their stars. This conclusion is consistent with the drop 
in the neutral baryon content of damped Lya galaxies at 2 > 3 (Storrie-Lombardi & Wolfe 1996), 
with the sharp decrease in quasar space density at z > 3 (cf. Kennefick et al. 1995), and with the 
morphology of z > 3 galaxies identified from their Lyman continuum break (Steidel et al. 1996a, b; 
Giavalisco et al. 1996). It would be extremely important to confirm the above results using a 
much larger statistical sample. 

4. The relative abundance patterns of the elements being studied here closely resemble those 
seen in Galactic halo stars. This provides (in our opinion) unambiguous evidence that the bulk of 
heavy elements in these high-redshift galaxies were produced by Type II supernovae, with little 
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contribution from stellar mass loss of low-to-intermediate mass stars or from Type la supernovae. 

The evidence for Type II supernova enrichment includes the overabundance of O relative to N, the 
overabundance of a-elements (Si, S) relative to the iron peak elements (Cr, Mn, Fe, Ni), and the 
underabundance of Mn relative to Fe and Al relative to Si (i.e., the odd-even effect). 

5. In agreement with earlier studies, we find Zn to be overabundant relative to Cr or Fe by 
an average of ~ 0.4 dex compared to their solar values. Although earlier studies have attributed 
this effect to selective depletion of Cr and Fe by dust grains, such an interpretation is inconsistent 
with many of the other elemental abundance ratios seen in these galaxies, most notably N/0 and 
Mn/Fe. We find no dependence of the measured elemental abundances on their condensation 
temperature in a way that would suggest the depletion by dust grains (cf. Jenkins 1987). The lack 
of significant variations in the S /Ni ratio across a large velocity interval (several hundreds of km 
s~^) in one damped Lya system also argues against the presence of dust in this galaxy. Therefore, 
we suggest that the overabundance of Zn relative to Cr in damped Lya galaxies may be intrinsic 
to their stellar nucleosynthesis. If this interpretation is correct, it will provide important new 
information to the theory of stellar nucleosynthesis. 

6. We discuss the physical conditions in damped Lya galaxies. The upper limits on the 
electron density in the absorbing gas as inferred from the C II*/C II ratios are consistent with 
the value found in the Galactic ISM. It is found that the absorption profiles of Al III in damped 
Lya galaxies always resemble those of the low-ionization lines. The profiles of Si IV and C IV 
absorption, while resembling each other in general, are almost always different from those of the 
low-ionization absorption lines. These results suggest that Al III is probably produced in the same 
physical region as the low-ionization species in the absorbing galaxies, while the high-ionization 
species (Si IV and C IV) mostly likely come from different physical regions. 

7. We discuss possible ways to get information on the history of star formation (i.e., 
continuous or episodic) in damped Lya galaxies, and on the shape of their stellar initial mass 
functions. 

8. We review the evidence for and against the hypothesis that damped Lya galaxies are disks 
or proto-disks at high redshifts (Wolfe et al. 1986), and discuss the implications. 

9. We determine upper limits on the temperature of the cosmic microwave background 
radiation at several redshifts using absorption from the fine structure level of the C II ion (i.e., C 
II*). These upper limits are consistent with the predicted increase of Tcmb with redshift. 
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Table 1. JOURNAL OF OBSERVATIONS 



Object 


em 


Date 


Coverage 


Exposure 


Q 0000-2620 


4.108 


13 Nov 95 


5107-7608 


6000s 






13 Nov 95 


5146-7660 


6000s 


Q 0216+0803 


2.992 


11 Oct 94 


4203-6633 


6000s 






27 Dec 94 


4223-6662 


6000s 






22 Aug 95 


5153-7670 


6600s 


Q 0449-1325 


3.097 


26 Dec 94 


4203-6632 


8000s 






26 Dec 94 


4224-6662 


8000s 


Q 0528-2505 


2.779 


27 Dec 94 


4541-7044 


6000s 






27 Dec 94 


6274-8697 


9000s 






23 Feb 95 


4540-7042 


6000s 






23 Feb 95 


6386-8695 


3000s 


Q 1425+6039 


3.173 


18 May 95 


4077-6540 


6000s 






18 May 95 


4059-6513 


8400s 






18 May 95 


3719-6179 


2400s 






19 May 95 


3719-6179 


9000s 






19 May 95 


3732-6199 


9000s 






19 May 95 


4077-6540 


2400s 


Q 1946+7658 


2.994 


12 Jun 94 


4013-6514 


4700s 






12 Jun 94 


4031-6541 


4000s 


Q 2212-1626 


3.992 


14 Nov 93 


5107-7608 


9000s 






15 Nov 93 


5146-7660 


9000s 






22 Aug 95 


4468-6887 


6000s 






22 Aug 95 


4492-6920 


6000s 






22 Aug 95 


5341-7775 


3000s 


Q 2231-0015 


3.018 


16 Jun 94 


4060-6515 


7600s 






17 Jun 94 


4032-6425 


6000s 


Q 2237-0608 


4.559 


10 Oct 94 


6401-8713 


12000s 






11 Oct 94 


6473-8801 


18000s 






20 Aug 95 


4899-7455 


12000s 






20 Aug 95 


4933-7350 


3000s 






21 Aug 95 


4933-7350 


9000s 



^ Redshift of the quasar as estimated from the Lya emission 
hne in the Keck spectrum except for Q 0528—2505, which is 
from Hewitt & Burbidge (1993). 
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Table 2. Measurements for the Zdamp = 3.3901 System Toward Q 0000-2620 



Ion 


A" 










Na{v) + a'^ 


adopt 


[Z/H]/ 


H I 


1215.670 


0.4164 










21.41+0.08 




Al II 


1670.787 


1.833 


-60 


+100 


0.286+0.020 


> 13.15 


> 13.15 


> -2.74 


Si II 


1526.707 


0.11 


-60 


+ 100 


0.327+0.008 


> 14.56 


> 14.56 


> -2.40 


Fe II 


1608.451 


0.0619 


-60 


+40 


0.178+0.013 


> 14.48 


> 14.48 


> -2.44 




1611.201 


0.00102 


-60 


+40 


0.003 + 0.017 


< 15.45 


< 15.45 


< -1.47 


C IV 


1548.195 


0.1908 


-170 


+130 


0.966+0.020 


> 14.73 


> 14.73 






1550.770 


0.09522 


-170 


+130 


0.814+0.0225 


> 14.869 







^Prom Morton 1991 unless revised values appear in the compilation of Tripp et al. 1996. 
^The velocity region over which the equivalent width and column density are measured. 
'^Rest-frame equivalent width and associated error. 

'^Column density from integrating the Na{v) profile. Upper limits are 4a". 
Tinal adopted column density. 

^Abundance relative to the solar value of Anders & Grevesse 1989. 

^The measurement is somewhat contaminated due to blending with the Mg II A2796 absorption 
at Zabs = 1-433. 
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Table 3. Measurements for the Zdamp = 2.2931 System Toward Q 0216+0803 



Ion 


A" 


r 




< 




Na{v) ± a'^ 


adopt 


[Z/H]/ 


H I 


1215.670 


0.4164 










20.45 + 0.16 




Si II 


1526.707 


0.11 


-340 


+140 


0.920+0.034 


> 15.06 


15.45+0.04 


-0.55 + 0.16 




1808.013 


0.00218 


-40 


+100 


0.153+0.012 


15.45 ± 0.04 






Zn II 


2026.136 


0.489 


-40 


+100 


0.067+0.029 


< 12.81 


< 12.81 


< -0.29 




2062.664 


0.256 


-40 


+100 


0.036+0.031 


< 13.11 






Pt- tt 
i^r 11 




n 1 

U. iUO 


/in 


^1 nn 

-T iUU 


U.UD4ztU.UZ4 


^ io.oy 


<^ lo.oy 


^ n 7/1 




2062.234 


0.078 


-40 


+100 


0.055+0.030 


< 13.61 






Fe II 


1608.451 


0.0619 


-80 


+140 


0.482+0.022 


> 14.81 


14.89+0.08 


-1.07 + 0.18 




2249.877 


0.00182 


-40 


+100 


0.027+0.028 


< 15.14 








2260.781 


0.00244 


-40 


+100 


0.067+0.026 


< 14.97 






Ni II 


1709.600 


0.0689 


-40 


+100 


0.064+0.024 


< 13.73 


13.60+0.07 


-1.10 + 0.17 




1741.549 


0.1035 


-40 


+100 


0.101+0.016 


13.60 ± 0.07 








1751.910 


0.06375 


... 9 


... 9 


... 9 


... 9 






Al II 


1670.787 


1.833 


-340 


+140 


1.091+0.030 


> 13.88 


> 13.88 


> -1.05 


AI III 


1854.716 


0.5602 


-340 


+140 


0.580+0.052 


13.75 ± 0.03 


13.74 ± 0.02 






1862.790 


0.2789 


-340 


+140 


0.359+0.033 


13.73 ± 0.02 






C IV 


1548.195 


0.1908 


-380 


+130 


1.460+0.032 


> 14.93 


> 15.06 






1550.770 


0.09522 


-380 


+130 


1.150+0.040 


> 15.06 







'^From Morton 1991 unless revised values appear in the compilation of Tripp et al. 1996. 
''The velocity region over which the equivalent width and column density are measured. 
^Rest-frame equivalent width and associated error. 

'^Column density from integrating the Na{v) profile. Upper limits are 4(7. 
^Final adopted column density. 

^Abundance relative to the solar value of Anders Sz Grevesse 1989. The uncertainty includes that from 
the Ar(H I). 

SBlended with the C IV A1550 absorption at Zabs = 2.7205. 
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Table 4. Measurements for the Zdamp = 1-7688 System Toward Q 0216+0803 



Ion 












Naiv) ± 


adopt 


[Z/H]/ 


H I 


1215.670 


0.4164 










20.00 + 0.18 




Si II 


1808.013 


0.00218 


-30 


+40 


0.044+0.011 


14.89+0.11 


14.89+0.11 


-0.66 ± 0.21 


Ni II 


1751.910 


0.06375 


-30 


+40 


0.020+0.005 


13.09+0.10 


13.09+0.10 


-1.16 + 0.21 


Mn II 


2576.877 


0.3508 


-30 


+40 


0.059+0.015 


12.48+0.11 


12.48+0.11 


-1.05 + 0.21 




2594.499 


0.2710 


-30 


+40 


0.052+0.017 


< 12.63 






Zn II 


2026.136 


0.489 


-30 


+40 


0.030+0.010 


< 12.34 


< 12.34 


< -0.31 




2062.664 


0.256 


... 9 


... 9 


... 9 


... 9 






Cr II 


2056.254 


0.105 


-30 


+40 


0.012+0.011 


< 13.06 


< 13.06 


< -0.62 




2062.234 


0.078 


-30 


+40 


0.033+0.014 


< 13.29 








2066.161 


0.0515 


-30 


+40 


0.011+0.020 


< 13.62 






Fe II 


2249.877 


0.00182 


-30 


+40 


0.013+0.010 


< 14.68 


14.53+0.09 


-0.98 ± 0.20 




2260.781 


0.00244 


-30 


+40 


0.035+0.007 


14.53+0.09 








2344.214 


0.1097 


-110 


+40 


0.395+0.013 


> 14.34 








2374.461 


0.0326 


-40 


+40 


0.253+0.011 


> 14.45 








2382.765 


0.3006 


-110 


+40 


0.533+0.011 


> 14.02 








2586.650 


0.0684 


-110 


+40 


0.407+0.024 


> 14.46 








2600.173 


0.2239 


-110 


+40 


0.534+0.028 


> 14.03 






Al III 


1854.716 


0.5602 


-150 


+40 


0.225+0.043 


13.23+0.08 


13.20+0.07 






1862.790 


0.2789 


-150 


+40 


0.102+0.022 


13.17+0.12 







^Prom Morton 1991 unless revised values appear in the compilation of Tripp et al. 1996. 
^The velocity region over which the equivalent width and column density are measured. 
^Rest-frame equivalent width and associated error. 

'^Column density from integrating the Na{v) profile. Upper limits are 4a". 
^Tinal adopted column density. 

^Abundance relative to the solar value of Anders & Grevesse 1989. The uncertainty includes that from 
the iV(H I). 

^Blended with an unidentified absorption feature. 



-48- 



Table 5. Measurements for the Zdamp = 1-2667 System Toward Q 0449-1326 



Ion 


A" 










Na{v) ± a'i 




Fe II 


2249.877 


0.00182 


-20 


+40 


0.082+0.005 


15.11+0.03 


15.13+0.02 




2260.781 


0.00244 


-20 


+40 


0.103+0.009 


15.14+0.03 






2344.214 


0.1097 


-60 


+50 


0.420+0.011 


> 14.27 






2374.461 


0.0326 


-40 


+50 


0.299+0.010 


> 14.67 






2382.765 


0.3006 


-60 


+140 


0.597+0.019 


> 13.98 






2586.650 


0.0684 


-60 


+50 


0.388+0.011 


> 14.39 






2600.173 


0.2239 


-60 


+50 


0.606+0.015 


> 14.02 




Mn II 


2576.877 


0.3508 


-20 


+40 


0.112+0.008 


12.86+0.03 


12.87+0.03 




2594.499 


0.2710 


-20 


+40 


0.097+0.015 


12.91+0.06 






2606.462 


0.1927 


-20 


+40 


0.099+0.021 


13.07+0.09 




Mg II 


2796.352 


0.6123 


-110 


+140 


1.080+0.015 


> 13.86 


> 14.11 




2803.531 


0.3054 


-110 


+140 


0.979+0.015 


> 14.11 




Mg I 


2852.964 


1.830 


-40 


+40 


0.267+0.009 


> 12.62 


> 12.62 



^Prom Morton 1991 unless revised values appear in the compilation of Tripp et al. 1996. 
^The velocity region over which the equivalent width and column density are measured. 
'^Rest-frame equivalent width and associated error. 
^Column density from integrating the Na{v) profile. 
^Final adopted column density. 

^No absolute abundance estimates are possible owing to the lack of A^(H I) information. 



Table 6. Measurements for the Zdamp = 1.1743 System Toward Q 0450-1312 



Ion 


A" 


J 




+ 




Nniv) ±a'^ 


adopt I'-'l ^'-i 


Mg II 


2796.352 


0.6123 


-50 


+260 


1.840 ± 0.038 


> 14.10 


> 14.10 




2803.531 


0.3054 


... 9 


... a 


... 9 


... 9 




Cr II 


2056.254 


0.105 


-20 


+130 


0.086 ± 0.022 


13.44 ± 0.12 


13.44 + 0.12 




2062.234 


0.078 


-20 


+130 


0.103 ± 0.032 


< 13.64 






2066.161 


0.0515 


-20 


+130 


0.044 ± 0.033 


< 13.83 




Mn II 


2576.877 


0.3508 


-20 


+130 


0.159 + 0.019 


12.96 + 0.08 


12.92 + 0.07 




2594.499 


0.2710 


-20 


+130 


0.100 + 0.020 


12.84 + 0.11 






2606.462 


0.1927 


-20 


+130 


0.069 ± 0.019 


< 13.02 




Fe II 


2249.877 


0.00182 


-20 


+130 


0.105 + 0.018 


15.06 + 0.10 


15.12 + 0.07 




2260.781 


0.00244 


-20 


+130 


0.129 ± 0.023 


15.17 + 0.09 






2344.214 


0.1097 


-50 


+240 


0.785 ± 0.042 


> 14.55 






2374.461 


0.0326 


-50 


+240 


0.494 ± 0.050 


> 14.86 






2382.765 


0.2006 


-50 


+240 


1.113 + 0.038 


> 14.29 






2600.173 


0.2239 


-50 


+240 


1.114 + 0.044 


> 14.33 




Zn II 


2026.136 


0.489 


-20 


+130 


0.067 ± 0.028 


< 12.80 


< 12.80 




2062.664 


0.256 


-20 


+130 


0.043 ± 0.035 


< 13.16 





^Prom Morton 1991 unless revised values appear in the compilation of Tripp et al. 1996. 
^The velocity region over which the equivalent width and column density are measured. 
^Rest-frame equivalent width and associated error. 

'^Column density from integrating the Na{v) profile. Upper limits are 4a". 
^Final adopted column density. 

^No absolute abundance estimates are possible owing to the lack of Ar(H I) information. 
SQnly partially covered. 
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Table 7. Measurements for the Zdamp = 0.8598 System Toward Q 0454+0356 



Ion 


A" 


r 




< 




Na{v) ± a'^ 


adopt 


[Z/H]/ 


H I 


1215.67 


0.4164 










20.76 ± 0.03 




Mg II 


2803.531 


0.3054 


-110 


+100 


1.450 + 0.017 


> 14.38 


> 14.38 


> -1.97 


Mg I 


2852.964 


1.830 


-90 


+60 


0.303 ± 0.027 


12.44 ± 0.04 


12.44 ± 0.04 




Cr II 


2056.254 


0.105 


-90 


+40 


0.093 + 0.108 


< 14.04 


< 14.04 


< -0.40 




2062.234 


0.078 


-90 


+40 


0.070 ± 0.098 


< 14.13 








2066.161 


0.0515 


-90 


+40 


0.043 ± 0.106 


< 14.34 






Mn II 


2576.877 


0.3508 


-90 


+40 


0.158 + 0.031 


12.96 + 0.09 


12.93 ± 0.06 


-1.36 + 0.07 




2594.499 


0.2710 


-90 


+40 


0.113 + 0.025 


12.90 + 0.09 








2606.462 


0.1927 


-90 


+40 


0.079 ± 0.033 


< 13.06 






Fe II 


2249.877 


0.00182 


-90 


+40 


0.074 ± 0.044 


< 15.33 


15.27 + 0.09 


-1.00 ± 0.09 




2260.781 


0.00244 


-90 


+40 


0.170 ± 0.034 


15.27 ± 0.09 








2344.214 


0.1097 


-100 


+80 


0.977 + 0.018 


> 14.74 








2374.461 


0.0326 


-100 


+80 


0.718 ± 0.024 


> 15.00 








2382.765 


0.2006 


-100 


+80 


1.141 + 0.016 


> 14.34 








2586.650 


0.0684 


-100 


+80 


1.014 + 0.016 


> 14.83 








2600.173 


0.2239 


-100 


+80 


1.228 + 0.011 


> 14.49 






Zn II 


2062.664 


0.256 


-90 


+40 


0.096 ± 0.095 


< 13.60 


< 13.60 


< +0.19 



^Prom Morton 1991 unless revised values appear in the compilation of Tripp et al. 1996. 
^The velocity region over which the equivalent width and column density are measured. 
'^Rest-frame equivalent width and associated error. 

'^Cohimn density from integrating the Na{v) profile. Upper limits are 4a". 
^Final adopted column density. 

^Abundance relative to the solar value of Anders 8z Grevesse 1989. 
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Table 8. Measurements for the Zdamp = 2.8110 System Toward Q 0528-2505 



Ion 


A" 




v_ 


v\ 




Na{v) + a'^ 


adopt 


[Z/H]/ 


n T 

J. J. J. 


l^ltJ.U 1 u 


41 f\A 










91 9f| _|_ n 1 fl 
ZjL.^yj u.xu 




C TT 


±00^.00^ 


1 978 


— 9fiO 


-l-'^^O 


9 904 -1- n 01 4S 

Zj.ZjO^ ^ U.UIt: 


^ 1 c: cog 


^ 1 "^8 


^ —9 1 R 
^ Zi. 10 


C TV 

KJ 1 V 


1 548 1Q5 


1Q08 


h 


h 


h 


h 


"> 1 5 00 






1 770 
±ooyj. I 1 u 


00^99 

U.Uc/OZZp 


— 180 

lOU 


-1-^90 


U.OO/ ZLL U.UOO 


~^ 1 l^ 00 
^ 10. UU 






N T 


1 1 QQ ^^n 
1 9nn 71 n 


1 *^97 

U. lOZf / 

u.uooo 
n 0449 


i 
i 

i 


i 
i 
i 


i 
i 

i 


i 
i 

i 






N V 

i>l V 


1 9*^8 S91 


1 ^7 

W. 1 (J 1 


— fiO 

\J\J 


-1-1 00 


1 fiq -1- 01 


1 4 nn -L n (-14 


1 9 qq _l_ f) fiA 






1 949 804 


0789^ 


— fiO 
uu 


-1-100 
n^iuu 


U.UOU ZLL U.Ul^ 


10.i/'4: ZIZ U.U/ 






O T 


1 ^09 1 fiQ 


04887 


— 9fi0 
^uu 


-1-^80 


1 QAc: -1- n 01 8J 


^ 1 ^ Q4J 
^ 10. 


^ 1 ^ Q4 


\ _9 in 


Mo- TT 


1 9'^Q Q9^ 


001 9^ 


—40 
^u 


-l-'^OO 
n^ouu 


n non _l n noo 

U.VJOVJ ZC U.VJOO 


\ 10.00 


\ 10.00 


<r —0 Q1 

\ U.t7l 




1 940 "^Q^ 


000fi9'i 


—40 


-l-'^OO 
n^ouu 


U.U^U ZC U.VJOO 


\ lU. liy 






Al TT 

jT-i 11 


1 fi70 7S7 
lU i u. / < 


x.ooo 


— 9fiO 
^uu 


-(-'^80 


1 7zt + 01 

Z(. 1 / 4: ZC VJ.VJIO 


^ 1 4 90 


^ 14 90 


^ _1 4S 


Al TTT 

Jr\.l 111 


1 8^4 71 fi 
100^. / lU 


'ifi09 


— 110 
1 lU 




U.c/UU zc u.uzu 


^ 1 4 04 


^ 1 4 07 






1 8fi9 7Q0 

lOU^. 1 ij\J 


978Q 


— 110 

1 IW 


-1-390 


n aAQ 4- fl 09^ 


^ 14 07 






cji TT 


1 9fi0 499 


1 007 
1 .uu 1 


— 9fi0 

ZrUU 


-1-^80 


9 97« -L fl 01 9^^ 


^ 1 4 70^^ 
^ 1^. / u 


lU. uu ZIZ U.U'4: 


_n 7"^ -1- 1 1 

U. f ZIZ U. 1 1 




1 ^04 ^70 


08fi 


— 1 90 


-l-'^80 


l.tJ^U ZIZ U.WIU 










1 ^9fi 707 


1 1 
u. 1 X 


— 9fi0 

^UU 


-l-'^80 


1 QQQ _L n noo/ 


^ 1 44^ 








1 808 01 


0091 8 


—80 
ou 


-l-'^OO 
n^ouu 


U.^t/O zc u.uz/ 


1 « nn _|_ n ri/i 
lu.uu zn u.u^ 






TV 

Ol 1 V 


1 ^Q^ 7^^ 


^14 


— 180 

lOU 


-1-^90 


1 nna -l n noo 

l.UUU ZIZ U.UZO 


^ 14 41 

^ 1^.4:1 


^ 14 "^9 






1 409 770 




— 180 

lOU 


-1-^90 


744 -1- 0^1 

U. i ZIZ U.UOl 


^ 14 ^9 






•=! TT 

O 11 


19^0 ^84 


00^4^^ 


—80 
ou 


-1-^00 


W.^Ul ZIZ \J.\JZ/ij 


10. uu in u.uo 


1 c: CQ _L n no 

10. Oc/ ZIZ U.UO 


_n 8^ -1- 10 

U.OO ZIZ U.IU 




1 9^"^ 81 1 

1 ZitJO . Oil 


01088 


—80 
ou 


-1-300 


f) 490 _|_ n nQfi 

yj.rt^ij _i_ u.uou 


1 c CO _|_ fl fiq 

-LtJ.t^O _1_ U.UO 








1 9^Q ^1 Q 


01 fi94 


m 


m 


m 


m 






Pr TT 

V^i 11 




10^ 


—40 
^u 


-l-'^OO 

T^OUU 


u.-LUc/ in u.u^u 


1 Q «C _1_ fl 1 9 


10. UO ZIZ U.lZi 


_-| 9^ -1- 1 

l.ZiO ZIZ U.IU 




zUdz.zo4 


U.U/o 


n 


n 


n 


n 








2066.161 


0.0515 


-40 


+300 


0.066 ± 0.037 


< 13.89 






Fe II 


1608.451 


0.0619 


-120 


+300 


1.264 + 0.018 


> 15.32 


15.45 + 0.11 


-1.26 + 0.15 




1611.201 


0.00102 


-40 


+300 


0.072 ± 0.026 


< 15.64 








2249.877 


0.00182 


-40 


+300 


0.194 + 0.091 


< 15.65 








2260.781 


0.00244 


-40 


+300 


0.274 + 0.066 


15.45 + 0.11 






Ni II 


1317.217 


0.146 


-40 


+300 


0.136 ± 0.025 


13.83 + 0.08 


13.89 + 0.03 


-1.56 + 0.10 




1370.132 


0.131 


-40 


+300 


0.159 ± 0.026 


13.92 + 0.07 








1454.842 


0.0596 


-40 


+300 


0.089 + 0.029 


< 14.02 








1709.600 


0.0689 


-40 


+300 


0.141 + 0.023 


13.94 + 0.07 








1741.549 


0.1035 


-40 


+300 


0.182 + 0.024 


13.87 + 0.06 
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Table 8 — Continued 



Ion 


A" 










Na{v) ± a"^ 


adopt 


[Z/H]/ 




1751.910 


0.06375 


-40 


+300 


0.124 + 0.025 


13.88 + 0.09 






Zn II 


2026.136 


0.489 


-40 


+300 


0.201 ± 0.030 


13.09 + 0.07 


13.09 + 0.07 


-0.76 + 0.12 




2062.664 


0.256 


o 


o 


o 


o 






Undetected Elements 


P II 


1532.533 


0.00761 


-80 


+300 


0.012 + 0.027 


< 14.84 


< 14.84 


< +0.07 


CI II 


1347.240 


0.118 


-80 


+300 - 


-0.012 ± 0.044 


< 13.97 


< 13.97 


< -0.50 


Ti II 


1298.697 


0.0952 


-80 


+300 


0.015 ± 0.031 


< 13.95 


< 13.95 


< -0.18 


Co II 


1466.203 


0.140 


-80 


+300 


0.022 + 0.031 


< 13.67 


< 13.67 


< -0.44 


Cu II 


1358.773 


0.380 


-80 


+300 


0.000 ± 0.030 


< 13.29 


< 13.29 


< -0.18 


Ga II 


1414.402 


1.801 


-80 


+300 


0.009 ± 0.032 


< 12.61 


< 12.61 


< +0.28 


Ge II 


1237.059 


0.876P 


-80 


+300 


0.015 + 0.040 


< 13.13 


< 13.13 


< +0.30 








For the absorption 


complext at v 


> 140 km 






C II 


1334.532 


0.1278 


+140 


+380 


0.904 ± 0.006 


> 15.18 


> 15.18 




C II* 


1335.708 


0.1149 


+140 


+300 


0.226 ± 0.010 


14.22 + 0.02 


14.22 ± 0.02 




N I 


1200.710 


0.0442 


+150« 


+300 


0.277 ± 0.009 


> 14.86 


< 15.33'' 




I 


1302.169 


0.04887 


+140 


+380 


0.807 + 0.007 


> 15.56 


> 15.56 




Si II 


1808.013 


0.00218 


+140 


+300 


0.129 + 0.012 


15.36 + 0.06 


15.36 + 0.06 




S II 


1250.584 


0.005453 


+140 


+300 


0.089 + 0.013 


15.10 + 0.07 


15.04 + 0.04 






1253.811 


0.01088 


+140 


+300 


0.129 + 0.015 


15.00 ± 0.05 






Fe II 


1608.451 


0.0619 


+140 


+300 


0.406 ± 0.008 


> 14.65 


> 14.65 






1611.201 


0.00102 


+140 


+300 


0.016 + 0.014 


< 15.37 








2260.781 


0.00244 


+140 


+300 


0.070 ± 0.036 


< 15.12 
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Note to Table 8: 

"Prom Morton 1991 unless revised values appear in the compilation of Tripp et al. 1996. 

''The velocity region over which the cqTiivalcnt width and column density are measured. 
'^Rest-frame equivalent width and associated error. 

'^Column density from integrating the Na{v) profile. Upper limits are 4cr. 
^Final adopted column density. 

■^Abundance relative to the solar value of Anders &; Grevesse 1989. The uncertainty includes 
that from the iV(H I). 

^May be slightly contaminated at the red edge by the C II* A1335 absorption. The effect on 
the column density limit should be negligible. 

'*Only partially covered by our spectrum. 

*The three N I lines are blended with each other. 

^Blended with the Si IV A1402 absorption at Zabs = 2.5381, but the effect on the column 
density limit should be negligible. 

^Contaminated at the blue edge by the S II A1259 absorption in the same system. The 
equivalent width is thus somewhat overestimated, but the effect on the column density limit 
should be negligible. 

'Slightly contaminated by the Al III A1854 absorption at Zabs = 2.1410, but the effects on the 
measured equivalent width and the column density limit should be negligible. 

"^Blended with the Si II A1260 absorption in the same system. 

"Blended with the Zn II A2062 absorption in the same system. 

"Blended with the Cr II A2062 absorption in the same system. 

^The /-value is from Cardelli et al. 1991a. 

^The f_ is adjusted to avoid the contaminating absorption at the left end. This should not 
have any significant effect on the estimated column density. 

''see text in S3. 5. 
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Table 9. Measurements for the Zdamp = 2.1410 System Toward Q 0528-2505 



Ion 


A" 












adopt 


Try /ttI f 

[Z/H]/ 


TT T 


l^lO.U / u 


n Al fizL 










on 7n 4- 08 
^u. 1 u m U.UO 




C TV 

± V 


1 ^48 1 

l*J^0.1c7«J 




1 AD 
— l^U 


-l-fiO 


n CI Q 4_ n 01 9 
u.uiy m u.uiz 


^ 14 l^7 


^ 1 4 71 
^ 1^. / 1 






1 ciKfi 770 

lOOU. / I u 




— 140 

— l^U 


4-fiO 

T^UU 


4QQ + 0014 


1 4 71 
^ 1^. / 1 






Al TT 


1 fi7n 787 
lU / u. / o / 


1 

1 .ooo 


— 1 90 

— izu 


-Ul 90 
izu 


n fil 9 -I- 01 ^ 
u.uiz zn U.UlO 


~^ 1 ^ 4fi 

^ lO.^U 


1 ^ 4fi 


^ — 1 79 


Al TTT 

X\L 111 


1 81^4 71 




^0 


n^ou 


1 04 4- 0079 
u.iu^ in u.uui 


■\o Qr. j_ r\ nog 
IZ.OO ZLL U.UO 


1 9 77 _|_ n n« 
iz. 1 / m U.UO 






1 8fi9 7Qn 


n 97RQ 


40 


n^ou 


04Q 4- OOQ 
u.u^c/ in u.uuy 


1 o 77 4- n 08 
1^. / 1 m U.UO 






Si TT 

Ol 11 


1 "^9^; 707 


1 1 


1 90 


4-1 90 


n KQc 4_ n ni 7 


^ 1 4 7Q 

^ 1^. ( a 


1 K 9« 4- n n/i 
lo.^u m u.u^ 


n QQ _|_ n OQ 




1 808 ni 

lOVJO.UlO 


n 0091 ^ 


fiO 


n^ou 


n nQ8 4- 008 
u.uyo m u.uuo 


1 c, 4_ fi n4 
lo.^u m u.u^ 






Cr TT 

KjL 11 






— ^0 


T^OU 


044 + 00'^ 
u.u^^ HI u.uuo 


1 Q flQ + 0'^ 
lo.ue? m U.UO 


1 Q 1 + 04 
lo.iu m u.u^ 


_i 9ft -U OQ 




OfifiO 9^4 


078 


— ^0 


T^OU 


n nor _|_ n nrir; 
u.uoo HI u.uuo 


1 Q 1 + Ofi 
lo.iu m u.uu 








9nfifi 1 fil 


0^1 ^ 


^0 


n^ou 


n 098 4- 007 
u.uzo in u.uui 


1 Q 1 « 4- n 1 1 

lo. lu m u. 1 1 






A/Tn TT 

IVlli 11 


9^7fi 877 


u.oouo 


■^0 
— ou 


4-'^0 
n^ou 


n 04^ 4- 01 1 
u.u^o in u.uii 


1 9 Qo 4- n 1 
i^.oo m u.iu 


1 9 QO 4- n 1 n 
i^.oo m u.iu 


i.oo zn u.±o 


Fe II 


1608.451 


0.0619 


-120 


+ 120 


0.320 ± 0.024 


> 14.68 


14.94 + 0.26 


-1.27 + 0.27 




1611.201 


0.00102 


-30 


+30 


0.017 + 0.010 


< 15.21 








2249.877 


0.00182 


h 


h 


h 


h 








2260.781 


0.00244 


h 


h 


h 


i 








2344.214 


0.1097 


-120 


+120 


0.778 + 0.019 


> 14.56 








2382.765 


0.2006 


-120 


+120 


1.091 + 0.024 


> 14.27 






Ni II 


1709.600 


0.0689 


-30 


+30 


0.019 + 0.009 


< 13.31 


13.22 + 0.06 


-1.73 + 0.10 




1741.549 


0.1035 


-30 


+30 


0.040 ± 0.007 


13.18 ± 0.08 








1751.910 


0.06375 


-30 


+30 


0.030 + 0.007 


13.27 + 0.10 






Zn II 


2026.136 


0.489 


i 


i 


i 


i 


< 12.28 


< -1.07 




2062.664 


0.256 


-30 


+30 


0.008 + 0.005 


< 12.28 







^Prom Morton 1991 unless revised values appear in the compilation of Tripp et al. 1996. 
^The velocity region over which the equivalent width and column density are measured. 
^Rest-frame equivalent width and associated error. 

^Column density from integrating the Na{v) profile. Upper limits are 4cr. 

*^Final adopted column density. 

^Abundance relative to the solar value of Anders k. Grevesse 1989. The uncertainty includes that from the 
A^(H I). 

^This line occurs in the red wing of the strong Si II A1526 line at Zahs = 2.8110 so it may be contaminated. 
^Blended with the Al III AA1854, 1962 lines at Zabs = 2.8110, respectively. 
'Blended with the Al II A1670 line at Zabs = 2.8110. 
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Table 10. Measurements for the Zdamp = 2.8268 System Toward Q 1425+6039 



Ion 


A" 








± a'' 




lyrc 
udopl 


[Z/II]/ 


H I 


1215.670 


0.4164 










20.30 ± 0.04 




C II 


1334.532 


0.1278 


-160 


+120 


0.842 ± 0.001 


> 15.19 


> 15.19 


> -1.67 


C IV 


1548.195 


0.1908 


-140 


+120 


0.545 ± 0.003 


14.32 ± 0.01 


14.32 ± 0.01 






1550.770 


0.09522 


-140 


+120 


0.322 ± 0.004 


14.31 + 0.01 






N I 


1199.550 


0.1327 


-120 


+60 


0.291 + 0.004 


> 14.62 


> 14.62 


> -1.73 




1200.223 


0.0885 


-120 


+60 


0.202 ± 0.004 


> 14.54 








1200.710 


0.0442 


... 9 


... 9 


... 9 


... 9 






Al II 


1670.787 


1.833 


-120 


+100 


0.634 ± 0.003 


> 13.55 


> 13.55 


> -1.23 


Si II 


1526.707 


0.11 


-50^^ 


+100 


> 0.471'^ 


> 14.78'* 


> 14.78 


> -1.07 


Si IV 


1393.755 


0.514 


-140 


+120 


0.288 ± 0.003 


13.62 + 0.01 


13.62 ± 0.01 






1402.770 


0.2553 


-140 


+120 


0.164 + 0.006 


13.62 + 0.01 






Fe II 


1608.451 


0.0619 


-120 


+100 


0.266 ± 0.004 


14.44 ± 0.01 


14.48 ± 0.04 


-1.33 + 0.06 




1611.201 


0.00102 


-40 


+40 


0.007 ± 0.002 


< 14.52 






Ni II 


1370.132 


0.131 


-40 


+40 


0.017 ± 0.001 


12.90 ± 0.03 


12.90 ± 0.03 


-1.65 + 0.00 




1454.842 


0.0596 


-40 


+40 


0.009 ± 0.002 


12.91 + 0.07 







^Prom Morton 1991 unless revised values appear in the compilation of Tripp et al. 1996. 
^The velocity region over which the equivalent width and column density are measured. 
^Rest-frame equivalent width and associated error. 

"^Column density from integrating the Na{v) profile. Upper limits are 4c7. 
^Final adopted column density. 

^Abundance relative to the solar value of Anders & Grevesse 1989. The uncertainty includes that from 
the iV(H I). 

^Contaminated by Lya forest absorption. 

'^The component at f ~ —90 km s^"*^ is blended with the C IV A1548 absorption at Zats = 2.7727. Thus 
we restrict the integration to v > —50 km s~^ to yield lower limits of the equivalent width and column 
density. 
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Table 11. Measurements for the Zdamp = 2.8443 System Toward Q 1946+7658 



Ion 


A" 






< 




Na{v) ± a'^ 


adopt 


[Z/H]/ 


H I 


1215.670 


0.4164 










20.27 ± 0.06 




C II 


1334.532 


0.1278 


-40 


+20 


0.119 ± 0.003 


> 14.14 


> 14.14 


> -2.69 


C II* 


1335.708 


0.1149 


-20 


+20 


0.001 ± 0.003 


< 12.76 


< 12.76 




C IV 


1548.195 


0.1908 


-100 


+110 


0.482 ± 0.006 


14.36 + 0.01 


14.36 ± 0.01 






1550.770 


0.09522 


-100 


+110 


0.319 + 0.006 


14.37 + 0.01 






N I 


1199.550 


0.1327 


-20 


+20 


0.011+0.003 


< 12.82 


< 12.82 


< -3.50 




1200.223 


0.0885 


-20 


+20 


-0.002 ± 0.003 


< 13.01 








1200.710 


0.0442 


-20 


+20 


0.008 ± 0.003 


< 13.28 






I 


1302.169 


0.04887 


-20 


+20 


0.081 ± 0.002 


> 14.46 


> 14.46 


> -2.74 


Al II 


1670.787 


1.833 


-40 


+20 


0.068 + 0.005 


12.26 + 0.03 


12.26 + 0.03 


-2.49 + 0.07 


Si II 


1304.370 


0.086 


-20 


+20 


0.040 ± 0.002 


13.64 ± 0.02 


13.63 ± 0.02 


-2.19 + 0.06 




1526.707 


0.11 


-20 


+20 


0.061 +0.002 


13.61 + 0.02 






Si IV 


1393.755 


0.514 


-100 


+110 


0.374 ± 0.008 


13.86 + 0.01 


13.86 + 0.01 






1402.770 


0.2553 


-100 


+110 


0.241 ± 0.007 


13.87 + 0.01 






S II 


1250.584 


0.005453 


-20 


+20 


0.003 ± 0.002 


< 14.00 


< 14.00 


< -1.54 


Fe II 


1608.451 


0.0619 


-20 


+20 


0.029 + 0.003 


13.38 + 0.04 


13.38 ± 0.04 


-2.40 ± 0.07 


Ni II 


1370.132 


0.131 


-20 


+20 


0.003 ± 0.002 


< 12.57 


< 12.57 


< -1.95 



"^From Morton 1991 unless revised values appear in the compilation of Tripp et al. 1996. 
^The velocity region over which the equivalent width and column density are measured. 
^Rest-frame equivalent width and associated error. 

'^Column density from integrating the Na{v) profile. Upper limits are Aa. 

*^Final adopted column density. 

^Abundance relative to the solar value of Anders Sz Grevesse 1989. The uncertainty includes that from 
the iV(H I). 
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Table 12. Measurements for the Zdamp = 1.7382 System Toward Q 1946+7658 



Ion 


A" 


r 






Wr ± (j"^ 


Na{v) ± a'^ 




Al III 


1854.716 


0.5602 


-20 


+40 


0.053 ± 0.004 


12.58 + 0.03 


12.58 + 0.03 




1862.790 


0.2789 


-20 


+40 


0.039 ± 0.004 


12.70 ± 0.04 




Si II 


1808.013 


0.00218 


-20 


+30 


0.033 ± 0.003 


14.76 ± 0.03 


14.76 + 0.03 


Cr II 


2056.254 


0.105 


-20 


+20 


0.025 ± 0.004 


12.84 + 0.07 


12.78 + 0.06 




2062.234 


0.078 


-20 


+20 


0.015 ± 0.003 


12.72 ± 0.08 






2066.161 


0.0515 


... 9 


... 9 


... 9 


... 9 




Fe II 


2249.877 


0.00182 


-20 


+20 


0.023 ± 0.003 


14.49 ± 0.05 


14.46 + 0.04 




2260.781 


0.00244 


-20 


+20 


0.026 ± 0.005 


14.41 ± 0.07 






2344.214 


0.1097 


-30 


+40 


0.254 ± 0.006 


> 14.17 




Zn II 


2026.136 


0.489 


-20 


+20 


0.006 ± 0.003 


< 11.83 


< 11.83 




2062.664 


0.256 


-20 


+20 


0.004 ± 0.003 


< 12.10 





"^From Morton 1991 unless revised values appear in the compilation of Tripp et al. 1996. 
^The velocity region over which the equivalent width and column density are measured. 
^Rest-frame equivalent width and associated error. 

'^Column density from integrating the Na{v) profile. Upper limits are Aa. 

'^Final adopted column density. 

^Abundance relative to the solar value of Anders & Grevesse 1989. The uncertainty includes that 
from the A^(H I). 

^Blended with the C IV A1548 absorption at Zabs = 2.6541. 
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Table 13. Measurements for the z^amp = 3.6617 System Toward Q 2212—1626 



Ion 


A" 










Na{v) + a'' 


adopt 


[Z/H]/ 


H I 


1215.670 


0.4164 










20.20 + 0.08 




C II 


1334.532 


0.1278 


-50 


+50 


0.276 ± 0.005 


> 14.53 


> 14.53 


> -2.29 


C II* 


1335.708 


0.1149 


-50 


+50 


0.017 + 0.007 


< 13.20 


< 13.20 




C IV 


1548.195 


0.1908 


-40 


+50 


0.140 ± 0.018 


13.73 + 0.06 


13.78 + 0.06 






1550.770 


0.09522 


-40 


+50 


0.101 + 0.016 


13.84 + 0.07 






I 


1302.169 


0.04887 


-50 


+50 


0.242 ± 0.001 


> 14.82 


> 14.82 


> -2.37 


Si II 


1304.370 


0.086 


-50 


+50 


0.087 ± 0.003 


13.91 + 0.01 


13.91 + 0.01 


-1.90 + 0.08 




1526.707 


0.11 


-50 


+50 


0.139 + 0.013 


13.92 + 0.04 






Si IV 


1393.755 


0.514 


... a 


... a 


... a 


... 9 


13.56 + 0.04 






1402.770 


0.2553 


-40 


+50 


0.111 + 0.009 


13.56 + 0.04 






Fe II 


1608.451 


0.0619 


-50 


+50 


0.038 + 0.034 


< 13.99 


< 13.99 


< -1.78 


Ni II 


1370.132 


0.131 


-50 


+50 


-0.012 + 0.012 


< 12.48 


< 12.48 


< -2.01 



^Prom Morton 1991 unless revised values appear in the compilation of Tripp et al. 1996. 
^The velocity region over which the equivalent width and column density are measured. 

'^Rest-frame equivalent width and associated error. 

"^Column density from integrating the Na{v) profile. Upper limits are 4(7. 
^Final adopted column density. 

^Abundance relative to the solar value of Anders & Grevesse 1989. The uncertainty includes that from the 
iV(H I). 

^Blended with the Mg II A2796 absorption line at Zabs = 1.3234. 
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Table 14. Measurements for the Zdamp = 2.0662 System Toward Q 2231-0015 



Ion 


A" 










Na{v) + a'^ 


adopt 


[Z/H]/ 


TJ T 
xl 1 


izio.D / 


n /1 1 RA 










on 1 n 1 n 




Al II 


1670.787 


1.833 


... 9 


... 9 


... 9 


... 9 






Al III 


1854.716 


0.5602 


-160 


+30 


0.173 + 0.012 


13.07 + 0.03 


13.07 ± 0.03 






1862.790 


0.2789 


-160 


+30 


0.084 + 0.017 


13.03 ± 0.09 






Si II 


1808.013 


0.00218 


-110 


+20 


0.091 + 0.012 


15.23 + 0.06 


15.23 + 0.06 


-0.88 + 0.12 


Cr II 


2056.254 


0.105 


h 


h 


h 


h 


< 13.26 


< -0.98 




2062.234 


0.078 


-110 


+20 


0.028 + 0.013 


< 13.26 








2066.161 


0.0515 


-110 


+20 


0.021 + 0.012 


< 13.40 






Fe II 


1608.451 


0.0619 


-160 


+30 


0.422 + 0.006 


> 14.72 


14.90 ± 0.18 


-1.17 + 0.20 




1611.201 


0.00102 


-110 


+20 


0.004 + 0.007 


< 15.08 






Ni II 


1709.600 


0.0689 


-110 


+20 


0.034 + 0.007 


13.30 + 0.10 


13.28 + 0.07 


-1.53 + 0.12 




1741.549 


0.1035 


-110 


+20 


0.048 + 0.009 


13.26 ± 0.08 








1751.910 


0.06375 


-110 


+20 


0.023 + 0.012 


< 13.45 






Zn II 


2026.136 


0.489 


-110 


+20 


0.035 + 0.008 


12.33 + 0.12 


12.33 + 0.12 


-0.88 ± 0.16 




2062.664 


0.256 


-110 


+20 


0.029 + 0.013 


< 12.74 







'^From Morton 1991 unless revised values appear in the compilation of Tripp et al. 1996. 
^The velocity region over which the equivalent width and column density are measured. 
^Rest-frame equivalent width and associated error. 

"^Column density from integrating the Na{v) profile. Upper limits are 4(7. 

*^Final adopted column density. 

^Abundance relative to the solar value of Anders & Grevesse 1989. The uncertainty includes that from 
the N(Kl). 

^Blended with the Si IV A1402 absorption at Zabs = 2.6523. 
'^Badly affected by a night sky line. 
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Tablc 15. Measurements for the z,i„„,p = 4.0803 System Toward Q 2237-0608 



Ion 


A" 






< 




Na{v) + a'' 




adopt 




[Z/H]/ 


H I 


1215. 


,670 


0.4164 










20, 


.52 + 0. 


,11 




C II 


1334. 


,532 


0.1278 


-140 


+80 


0.579 ± 0.002 


> 14.89 




> 14. 


,89 


> -2.19 


C II* 


1335. 


,708 


0.1149 


-130 


+60 


-0.001 ± 0.003 


< 12.83 










C IV 


1548. 


,195 


0.1908 


-120 


+80 


0.198 + 0.019 


13.81 + 0.04 


13, 


.81 + 0. 


,04 






1550. 


,770 


0.09522 


-120 


+80 


0.109 + 0.019 


13.81 + 0.08 










Al II 


1670. 


,787 


1.833 


-130 


+60 


0.243 ± 0.013 


12.85 + 0.02 


12 


.85 + 0. 


,02 


-2.15 + 0.11 


Si II 


1526. 


,707 


0.11 


-130 


+60 


0.291 + 0.007 


14.27 + 0.02 


14, 


.27 + 0. 


,02 


-1.80 + 0.11 


Si IV 


1393. 


,755 


0.514 


-120 


+80 


0.183 + 0.010 


13.41 + 0.03 


13, 


.44 + 0. 


,03 






1402. 


,770 


0.2553 


-120 


+80 


0.119 + 0.006 


13.47 + 0.03 










Fe II 


1608. 


,451 


0.0619 


-130 


+60 


0.087 + 0.022 


13.85 + 0.11 


13, 


.85 + 0. 


,11 


-2.18 + 0.15 


Ni II 


1370. 


,132 


0.131 


-100 


+300 


-0.003 + 0.008 


< 13.17 




< 13. 


,17 


< -1.6 



^Prom Morton 1991 unless revised values appear in the compilation of Tripp et al. 1996. 
^The velocity region over which the equivalent width and column density are measured. 

^Rest-frame equivalent width and associated error. 

"^Column density from integrating the Na{v) profile. Upper limits are Aa. 
^Final adopted column density. 

^Abundance relative to the solar value of Anders & Grevesse 1989. The uncertainty includes that from 
the Ar(H I). 
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Table 16. Dummy Table as Placeholder 



Object Zd^rnp N{C II) iV(C II*)-^ (cm'^) reference 



Table 17. Dummy Table as Placeholder 



Object Zdamp N{C II) N{C 11*)" Ue (cm-3) reference 
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Table 18. Electron Densities 



Object 


^damp 


N{C II) 


N{C II*)" 


Ue (cm ^) 


reference 


0528-2505 


2.8110^ 


> 15.18 


14.22 


< 1.6 


1 


1946+7658 


2.8443 


> 14.14 


< 12.46 


< 0.30 


1 


2212-1626 


3.6617 


> 14.53 


< 12.90 


< 0.34 


1 


2237-0608 


4.0803 


> 14.89 


< 12.53 


< 0.06 


1 


1202-0725 


4.3829 


> 14.96 


< 13.06 


< 0.18 


2 


References. 


— (1) this paper; (2) Lu et al. 


1996a. 







^ Upper limits are 2a. 



For the components at v > 140 km s 
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Table 19. Comparisons of Different Ionization Species 



(Jbject 


^damp 


A 1 TT / A 1 TTTa 
Ai 11/ Ai 111 


Low ion vs AI ill 


bi IV vs O IV 


Low Ion vs High Ion 


Kel 


nnnn 9fi9n 
uuuu — ZOZ\J 


o.oyui 








QmerenTj 


1 
1 


0100+1300 


2.3090 






similar? 


different 


2 


0201+3634 


2.4620 


» 1 




similar 


different 


3 


0216+0803 


2.2931 


» 1 


similar 




different? 


1 




1.7688 




similar 






1 


0528-2505 


2.8110 


» 1 


similar 


similar 


similar? 


1 




2.1410 


» 1 


similar 




different 


1 


1202-0725 


4.3829 






similar 


different 


4 


1331+1704 


1.7764 








different 


5 


1425+6039 


2.8268 






similar? 


different 


1 


1946+7658 


2.8443 






similar 


different 


1 




1.7382 




similar 






1 


2212-1626 


3.6617 






similar 


different 


1 


2231-0015 


2.0662 




similar 






1 


2237-0608 


4.0803 






similar? 


different 


1 



References. — (1) this paper; (2) Wolfe et al. 1994; (3) Prochaska & Wolfe 1996; (4) Lu et al. 1996a; (5) 
Wolfe 1995. 

'^This column compares the relative strength of the Al II A1670 absorption with that of the Al III AA1854, 
1862 absorption. A >> 1 means the Al II absorption is much stronger than the Al III absorption. 

'^Thcsc columns compare the absorption profiles of low ion lines with Al III, and with high ion lines (ic, Si 
IV and C IV). The entry "similar" means that the two ions being compared have similar absorption profiles 
in terms of the velocity distribution and relative strength of the components. The entry "similar?" means 
the two ions being compared are more similar than different. The entries "different" and "different?" should 
be interpreted in an analogous way. 
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Tablc 20. Measurements of Cosmic Microwave Background Temperature 



Object 


redshift 


Tel (K) 


TcMB (K) 


reference 


1331+1704 


1.7760 


10.4 + 0.5 


7.6 


1 






7.4 + 0.8 


7.6 


1 


0528-2505 


2.8110 


< 31.5^ 


10.4 


2 


1946+7658 


2.8443 


< 20.1 


10.5 


2 


0636+6801 


2.9090 


< 16.0'^ 


10.7 


3 


2212-1626 


3.6617 


< 20.6 


12.7 


2 


2237-0608 


4.0803 


< 14.9 


13.9 


2 


1202-0725 


4.3829 


< 18.1 


14.7 


4 



References. — (1) Songaila et al. 1994a; (2) this paper; (3) Songaila et al. 1994b; (4) Lu et 
al. 1996a. 

^ Excitation temperature as derived either from N{C U*)/N{C II) or from N{C I*)/N{C I). 
Upper hmits are 2a. 

^ Predicted temperature of the cosmic microwave background from 2.73(1+2;). 
For the absorption complext at v > 140 km (see Table 8). 

Songaila et al. (1994b) quoted T^x < 13.5 K using N{C II)=4.6 x 10^^ cm-^ estimated from 
profile fitting. Given the heavy saturation of the C II A1334 absorption, we consider this N{C 
II) very uncertain, and prefer to use N{C II)> 1.5 x 10^^ cm"^ estimated from the observed 
equivalent width of the C II absorption. 
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Fig. 1. — Damped Lya absorption in the Zdamp = 3.3901 system toward Q 0000—2620. The sohd 
curve represents the best-fit damping profile with iV(H I)= 2.6 x 10^^ cm~^. The dotted curves are 
damping profiles with iV(H I)= (2.2,3.1) x lO^^ cm'^. 

Fig. 2. — Velocity profiles of metal lines in the Zdamp = 3.3901 damped Lya system toward Q 
0000—2620. The zero velocity is fixed at z = 3.39005. Absorption features marked with "+" signs 
are unrelated to the absorption line in question. Lines that occur in the Lya forest or are blended 
with other absorption lines are indicated as such. 

Fig. 3. — Velocity profiles of metal lines in the Zdamp = 2.2931 damped Lya system toward Q 
0216+0803. The zero velocity is fixed at z = 2.29310. The absorption at v > 130 km in the 
C IV A1548 panel (marked with a "+" sign) is from C IV A1550, and that at v < —380 km s~^ in 
the C IV A1550 panel (marked with a "+" sign) is from C IV A1548. 

Fig. 4. — Velocity profiles of metal lines in the Zdamp = 1-7688 damped Lya system toward Q 
0216+0803. The zero velocity is fixed at z = 1.76880. 

Fig. 5. — Velocity profiles of metal lines in the Zdamp = 1.2667 damped Lya system toward Q 

0449- 1326. The zero velocity is fixed at z = 1.26666. The strongly saturated Fe II A2600 hne 
is not shown due to space considerations. Features marked with "+" signs are unrelated to the 
absorption line in question. 

Fig. 6. — Profiles of apparent column densities, Na{v), for the Mn II lines in the Zdamp — 

1.2667 

damped Lya system toward Q 0449—1326. 

Fig. 7. — Velocity profiles of metal lines in the Zdamp = 1-1743 damped Lya system toward Q 

0450- 1312. The zero velocity is fixed at ^ = 1.17428. The data are from Churchill & Vogt (1997). 

Fig. 8. — Velocity profiles of metal lines in the z^amp = 0.8598 damped Lya system toward Q 
0453+0356. The zero velocity is fixed at z = 0.85980. The data are from Churchill & Vogt (1997). 

Fig. 9. — Velocity profiles of metal lines in the Zdamp = 2.8110 damped Lya system toward Q 
0528—2505. The zero velocity is fixed at z = 2.81100. Features marked with "+" signs are 
unrelated to the absorption in question. Lines that are blended with other absorption lines are 
indicated as such. 

Fig. 10. — Velocity profiles of metal lines in the Zdamp = 2.1410 damped Lya system toward Q 
0528—2505. The zero velocity is fixed at z = 2.14100. Absorption features marked with "+" signs 
are unrelated to the absorption line in question. The Ni II A1709 and Cr II A2066 lines, though 
detected, are not shown due to space considerations. Lines that are blended with other absorption 
lines are indicated as such. 

Fig. 11. — Lya and Ly/3 absorption line profiles in the z^amp = 2.8268 damped Lya system toward 
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Q 1425+6039. The dotted curves show the damping profiles with Zabs = 2.8268, iV(H I)= 2 x 10^° 
cm~^, and Doppler b = 10 km s~^. Note the excess absorption at the red wing of the line profiles, 
which comes from the absorbing subcomplex near Zabs = 2.83058 {v ~ 300 km s~^ in figure 13). 

Fig. 12. — Same as figure 11, but with one additional component added to the fit at Zabs = 2.83058 
with Ar(H I)= 1 X 10^9 cm-2 and 6 = 10 km s'^ 

Fig. 13. — Velocity profiles of metal hues in the z^amp = 2.8268 damped Lja system toward Q 
1425+6039. The zero velocity is fixed at z = 2.82680. Absorption features marked with "+" signs 
are unrelated to the absorption line in question. Lines occurring in the Lya forest are indicated 
as such. Note the expanded vertical scales for some panels. Essentially all absorption displayed in 
the C 11* A1335 panel is from the C 11 A1334 line. Only the components at v < 120 km will be 
considered as part of the damped Lya system (see main text). 

Fig. 14. — Damped Lya absorption in the Zfiamp — 

2.8443 system toward Q 1946+7658. The soHd 
curve represents the best-fit damping profile with A^(H 1)=1.86 x 10^° cm~^. The dotted curves 
are damping profiles with A^(H 1)=(1.62, 2.14) x 10^° cm~^. 

Fig. 15. — Velocity profiles of metal lines in the z^amp = 2.8443 damped Lya system toward Q 
1946+7658. The zero velocity is fixed at z = 2.84430. Absorption features marked with "+" signs 
are unrelated to the absorption line in question. Lines occurring in the Lya forest are indicated as 
such. 

Fig. 16. — Velocity profiles of metal lines in the Zdamp = 1.7382 damped Lya system toward Q 
1946+7658. The zero velocity is fixed at ^; = 1.73820. Absorption features marked with "+" signs 
are unrelated to the absorption line in question. Lines blended with other lines are indicated as 
such. 

Fig. 17. — Damped Lya absorption in the Zdamp = 3.6617 system toward Q 2212—1626. The solid 
curve represents the best-fit damping profile with A^(H I)=1.6 x 10^° cm~^. The dotted curves are 
damping profiles with N{B I)=(1.3, 1.9) x 10^° cm'^. 

Fig. 18. — Velocity profiles of metal lines in the Zdamp = 3.6617 damped Lya system toward Q 
2212—1626. The zero velocity is fixed at z = 3.66170. Absorption features marked with "+" signs 
are unrelated to the absorption line in question. Lines blended with other absorption lines are 
indicated as such. 

Fig. 19. — Velocity profiles of metal lines in the z^amp = 2.0662 damped Lya system toward Q 
2231—0015. The zero velocity is fixed at z = 2.06615. Absorption features marked with "+" signs 
are unrelated to the absorption line in question. 

Fig. 20. — Damped Lya absorption in the Zdamp = 4.0803 system toward Q 2237—0608. The solid 
curve represent the best-fit damping profile with A^(H 1)= 3.0 x 10^° cm~^. The dotted curves are 
damping profiles with Ar(H I)=(2.3, 3.9)xl020 cm'^. 
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Fig. 21. — Velocity profiles of metal lines in the Zdamp = 4.0803 damped Lya system toward Q 
2237—0608. The zero velocity is fixed at z = 4.08026. Absorption features marked with "+" signs 
are unrelated to the absorption line in question. 

Fig. 22. — Age-metallicity relations. The conversion between age of the universe and redshift is for 
the cosmological parameters indicated in the figure. The solid circles are for damped Lya galaxies 
in Table 16. The "+"s are for Galactic disk stars in the sample of Edvardsson et al. (1993). 

Fig. 23. — Abundance ratios of selected elements against [Fe/H] for the sample of damped Lya 
galaxies given in Table 16. In a few cases, except for the [Cr/Fe] vs [Fe/H] panel, Cr abundance 
has been used in place of Fe abundance when the latter is not available. 

Fig. 24. — Similar to fig. 23, but for diffuse ISM clouds in the Galactic disk and halo as given in 
Table 17. Note that the spread in [Fe/H] in the horizontal axis reflects the different levels of Fe 
depletion in the ISM clouds rather than intrinsic differences in the clouds' metallicity. 

Fig. 25. — Similar to fig. 23, but for Milky Way disk and halo stars. These ratios reflect the 
chemical enrichment processes during the past history of the Milky Way. Sources of the data are 
given in §4.2. 

Fig. 26. — Elemental abundances as a function of condensation temperature for damped Lya 
galaxies. Note that the abundances are plotted relative to the Fe abundance (or Cr abundance, in 
case Fe measurement is not available) in the same system. For Al, Si, S, Cr, Mn, Fe, Ni, and Zn, 
the average values of the actual measurements (i.e., excluding upper and lower limits) are plotted 
(solid dots), with the vertical bars indicating the range of the measured values. For C, N, and O, 
individual limits are plotted since no actual measurements are available. Lower limits are indicated 
by open triangles, and upper limits are indicated by inverted solid triangles. 

Fig. 27. — The three left panels show the profiles of apparent column density per unit velocity, 
Na{v), for S II, Si II, and Ni II in the Zdamp = 2.8110 damped Lya system toward Q 0528—2505, 
based on the average of two S II lines, the Si II A1808 line, and the average of six Ni II lines. The 
three panels to the right give the Na{v) ratios for S Il/Si II, S II/Ni II, and Si II/Ni II. 

Fig. 28. — Temperature of the cosmic microwave background as inferred from observations of high- 
redshift quasar absorption systems. The dotted line is the predicted relation of Tqmb = 2.73(1 + 2;) 
from Big Bang cosmology. 
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TABLE 16 

Elemental Abundances of Damped Lyman-alpha Galaxies 
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Notes to Table 16: 

" Taking the S II column density discussed in §3.1 as an upper limit given the possibility of contamination by Lya forest absorption. 

' The #(II I) of the system is unknown (see discussion in §3). The elemental abundances are given for an assumed [Fe/II]= —1.50, which is typical of the 
sample galaxies. 



TABLE 17 



Elemental Abundances of Diffuse ISM Clouds 
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-0.07 


3 




B 


-0.12 




+0.19 




-0.39 


-0.20 




3 


HD 167756 


A+B+C 


> -0.53 




-1.21 




-1.43 




-0.22 


4 


HD 93521 


1 


-0.03 


+0.31 




-0.29 


-0.38 






5 




3 


-0.16 


+0.06 




-0.58 


-0.75 






5 




4 


-0.31 


-0.13 




-0.64 


-0.61 






5 




6 


-0.29 


-0.04 




-0.73 


-0.76 






5 







n zLQ 

— u .^y 


— u .uo 




1 1 1 
— 1.11 


1 90 
— 1 .zu 








72 Vel 


6 


-0.63 


+0.03 




-1.02 


-1.59 






6 




7 


-0.65 


+0.03 




-0.65 


-1.12 






6 


HD 149881 


2d 


-0.22 




-0.54 


-0.52 


-0.72 




-0.11 


7 




gd 


-0.47 




-0.88 


-0.88 


-0.93 




-0.11 


7 




gd 


-0.29 




-0.67 


-0.67 


-0.83 




-0.11 


7 


HD 116852 


Group 2« 


> -1.00 


-0.03 


-0.63 


-0.63 


-0.68 


-0.92 




8 




Group 3* 


> -0.49 


-0.13 


< -0.47 


-0.56 


-0.77 


-0.80 




8 


HD 18100 


[-60,60] 


-0.23 


+0.03 


-0.72 


-0.66 


-0.80 


-1.15 




8 


^ Per 




> -1.26 




-2.08 


-1.31 


-2.10 


-2.45 


-0.56 


9 



REFERENCES: (1) Savage, Cardelli, k Sofia 1992; (2) Sofia, Cardelli, k Savage 1994; (3) Sofia, Savage, 
k Cardelli 1993; (4) Cardelli, Sembach, k Savage (1994); (5) Spitzer k Fitzpatrick 1993; (6) Fitzpatrick k 
Spitzer 1994; (7) Spitzer k Fitzpatrick 1995; (8) Sembach k Savage 1996; (9) Cardelli et al. 1991b. 



Notes to Table 17: 

" This refers to the component name or number given in the original paper. In the case of ^ Per, the 
measurements refer to the entire absorption. 

' Reference 1 assigned a log#(H I)=19.74 to component A, while reference 2 used log#(H I)=20.16. The 
values given here are for log#(H I)=20.16. 

Reference 1 assigned a log #(H I)=21.12 to component B, while reference 2 used a log #(H I)=20.09. 
Here we adopt log#(H I)=20.09. 

Reference 7 found [Zn/H]=— 0.11 ± 0.08 for the entire absorption, and gave depletions of other elements 
relative to Zn in each component in order to avoid the difficulty of assigning #(H I) to each component. 
Here we assume [Zn/H]=— 0.11 for each individual components, and convert the depletions relative to Zn 
back to abundances relative to hydrogen. 

* Due to difficulties of assigning #(H I) to each group, the abundance measurements in reference 10 were 
given relative to Zn assuming Zn is undepleted in these warm clouds. 



